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ABSTRACT 

(Distribution  Limitation  Statement  A) 


This  study  presents  a  study  of  the  polarization  of  infrared  waves  due  to 
propagation  through  the  atmosphere.  Cumulus  cloud  and  Haze  M  particle 
size  distributions  are  considered  for  various  optical  depths  at  wave¬ 
lengths  of  4.0  and  10., 0y.  Results  at  these  wavelengths  are  compared  ,to 
results  at  a  wavelength  of  0.5y  for  the  same  physical  atmosphere. 
Initially  the  radiative  transfer  integral  equations  are  derived  which 
include  the  Stokes  polarization  vector.  From  these  integral  equations 
a  Monte  Carlo  simulation  technique  is  used  to  simulate  photon  histories 
through  the  atmosphere.  The  effects  of  multiple  scattering  on  the 
polarization  state  of  the  scattered  intensity  Is  included.  A  diffuse 
scattering  ground  which  causes  complete  depolarization  of  the  incident 
intensity  is  included  in  the  simulation  model.  Numerical  results  are 
obtained  for  cloudy  and  hazy  atmospheres.  It  is  shown  that  the  degree 
of  depolarization  is  a  function  of  atmospheric  particle  density  and: 
wavelength,  thereby  leading  to  the 'conclusion  that  in  a  linear*  homo¬ 
geneous,  .sotropic  scattering  medium,  depolarization  is  a  result  of 
multiple  scattering.  Significant  depolarization  at  wavelengths  of 
0.5  and  4:0y  is  shown  to  occur  for  cumulus  clouds  and  optical  depths  of 
3.0  when  the  cloud  is  1  km- in  height.  Polarization  factors  are  obtained 
with  and  without  ground  reflections,  with  ground  albedos  playing,  a 
significant  part  in  the  polarization  of  the  backs catter  intentity. 
Detectors  are  .judiciously  placed  .below  and  above  the  cloud  layer  so  that 
forward  and  backscatter  are  measured. 
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SECTION  I 


INTRODUCTION 

Polarization  is  a  fundamental  characteristic  of  electromagnetic 
waves.  The  polarization  of  a  wave  propagating  in  a,  media  is  defined 
as  the  orientation  and  direction  of  the  electric  field' vector.  The 
polarization  properties  of  electromagnetic  waves  have  been  used  spar¬ 
ingly  in,  describing  the  properties  of  a  scatterer.  In  radar  circular 
polarization  has  been  used  to  discriminate  between  desired  targets  and 
rain.  Orthogonal  polarization  has  been  used  in  communication  systems 
to  prevent  interference  between  two  overlapping  channels. 

One  of  the  obstacles  in  the  greater  use  of  polarization  has  been 
the  difficulty  of  controlling  and  measuring  the  polarization  state 
of  a  system.  Interference  between  targets  which  causes  depolariza¬ 
tion  beyond  that  which  was  expected  has  been  another  problem.  The 
polarization  of  television  signals  used  in  the  United  States  is 
different  from  that  used  in  Great  Britain.  In  Great  Britain  the 
polarization  is  vertical.  Vertical  polarization  is  used  to  reduce 
ground  reflections  and  interference  from  overhead  wires.  Horizontal 
polarization  is  used  in  the  United  States.  Ignition  and  industrial 
noise  contains  vertical  components  which  are  usually  greater  than  the 
horizontal  components  and  therefore  interfere  with  vertically  polarized 
signals. 

Optimum  discrimination  can  be  achieved  if  the  unwanted  reflected 
wave  is  completely  depolarized  with  respect  to  the  wave  reflected  from 
the  desired  target.  In  practice  this  situation  is  never  achieved.  In 
real  life  both  the  desired  and  unwanted  signals  will  be  partially 
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polarized  with  the  extent  of  depolarization  of  one  greater  than  the 
other.  If  the  difference  in  depolarization  between  the  two  signals  is 
sufficient  some  degree  o?  discrimination  can  be  achieved. 

Recently,  Manz  [l]1  has  proposed  a  discrimination  technique  using 
polarization  which  will'  separate  simultaneous  returns  from  cloudc.  and 
terrain  when  using  a  laser.  In  his  model  of  the  cloud  and  the  terrain, 
Manz  has  assumed  that  the  cloud  is  composed  of  Mie  particles,  i.e., 
spherical  particles,  and,  that  the  terrain  is  rough  for  all  wavelengths 
of  interest.  The  backscattered* return  from  the  cloud  consisting  of 
Mie  particles  will  be  depolarized  less  than  the  backsoattered  return 
from  the  terrain.  Mie  s  work  [2]  shows  that  for  incident  linear 
polarization,  as  Manz  has  proposed,  perpendicular  or  parallel  to  the 
-plane  of  scatter,  no  depolarization  will  occur.  On  the  other  hand, 
for  very  rough  terrain  one  would  expect,  on  the  average,  to  receive 
equal  parallel  and1  perpendicular  components  for  an  incident  linearly 
polarized  wave.  Recent  workers  [3,  4,  5]  have  calculated  the  cross- 
polarization  ratio  for  rough  terrain  defined  by 


where  Ec  is  the  cross  polarized  electric  field  component  of  the  scat¬ 
tered  electromagnetic  field  and  E  is  the  component  of  the  scattered 

r 

electromagnetic  field  which  is  colinear  with  the  incident  field.  The 
asterisk  denotes  the  complex  conjugate.  Using  random  sampling  from  a 
normal  distribution  to  <\iiulate  the  surface  function,  Beckm&nn  [6]  shows 


* 


^he  numbers  in  brackets  refer  to  the  corresponding  numbers  in  the 
references. 
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that  the  cross-polarization  ratio  Q  depends  on  the  angle  of  incidence* 
the  rms  slope  of  the  terrain,  and  the  dielectric  const'nt.  Minimum' 
depolarization,  according  to  Beckmann  and  others  [3,  4,  5],  occurs  at  a 
zero  angle  of  incidence.  The  ratio  Q  increases  with  angle  of  incidence. 
It  is  therefore  reasonable  to  assume  that  at  optical  frequencies,  where 
the  terrain  is  relatively  rougher  than  at  microwave  frequencies-,  the 
ratio  Q  is  equal  to  at  least  one-half.  While  ft  is  true  that,  at 
optical  frequencies,  the  distance  between  particles  in  a  cloud  is  large 
compared  to.  a  wavelength  so  that  near  field  effects  may  be  neglected, 
far-field  multiple  scattering  will  probably  cause  some  depolarization 
of  the  backscattered  wave.  Therefore  Manz's  discrimination  technique 
does  not  give  a  sure  way  of  separating  returns-  from  clouds  and  rough 
targets  for  all  situations  as  borne  out  by  his  experimental  results. 

Fung  has  also  shown  that,  for  smooth  reflecting  plane  surfaces;,  the 
cross-polarized  component  will  be  zero  for 

(a)  locally  grazing  incidence, 

,(b)  incident  polarization  parallel  with  or  perpendicular  to  the 
local  plane  of  incidence,  and 

(c)  R+  +  R  ■»  0.  Where  R+  is  the  Fresnel  reflection  coefficient 
for  incident  polarization  in  the  plane  of  incidence  and  R"  is 
the  Fresnel  reflection  coefficient  for  incident  polarization 
perpendicular  to  the  plane  of  incidence. 

Depolarization  of  electromagnetic  waves  can  be  caused  by  different 
phenomena.  As  implied  above,  a  wave  can  be  depolarized  by  multiple 
scattering.  One  has  to  take  care  when  discussing  multiple  scattering 
since  it  may  be  understood  to  occur  when  two  objects  are  within  mutual 
near-fields  cr  it  may  be  understood  to  occur  as  a  result  of  multiple 


reflections,  i.e.,  the  reflecting'  objects,  may  not  be  in  mutual  near¬ 
fields.  One  example  of  depolarization  due  to  near-field  scattering  is 
that  of  two  or  more  spheres  in  close  proximity,  where  the  individual 
scattering  functions  can  no  longer  be  considered  separately  but  must  be 
derived  only  by  including  the  interdependence  of  the  individual  sphere's 
field.  In  this  case  the  conglomerate  scattering  function  will  give 
rise,  in  general,  to  polarization  in  a  particular  direction  which  is 
different  than  if  one  had  considered  the  individual  scattering  functions 
separately.  When  considering  depolarization  due  to  multiple  reflection 
or  multiple  bounce  one  assumes  that  the  scatterers  are  in  each  other's 
radiation  or  far-field  zone.  In  this  case  one  can  consider  the  indi¬ 
vidual  scattering  functions  separately;  and  these  will,  in  general, 
lead  to  a  different  polarization  state  than  had  the  scatterers  been  in 
close  proximity.  The  difference  in  polarization  state- between  the  two 
types  of  multiple  scattering  may  be  greater  than  or  less  than  that 
given  by  the  other  method,  depending  on  the  scattering  angle  being  con¬ 
sidered.  There  is  one  possible  solution  to  this  dilemma.  If  we  approach 
the  problem  of  scattering  from  a  photon-particle  collision  standpoint 
and  if  we  further  select  the  particle-scattering  function  which  was 
determined  for  a  single  particle  in  the  radiation  zone  to  represent  the 
probability  density  for  scattering  in  a  certain  direction,  we  can  then 
bring  particles  as  close  together  as  we  wish  without  touching.  This 
approach  will  give  some  approximation  to  the  multiple  scattering  prob¬ 
lem  since  if  we  run  many  photon-particle  collision  histories  for  the 
conglomerate,  we  may  now  measure  the  probability  scattering  function 
for  the  conglomerate.  The  new  probability  scattering  function  will  not 
in  general  be  equal  to  the  sum  of  individual  particle  scattering  func¬ 
tions  because  of  the  multiple  collisions  which  were  allowed  to  take 
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place.  For  example*  cohside  the  scattering  of  an  electromagnetic 
wave  by  a  sphere.,  ^single  sphere  may  have  a  scattering  function  which 
is  highly  peaked  in  the  forward  scattering  direction  (figure  1). 


Figure  1.  Scattering  Function  for  a  Single  Sphere 

For  this  particular  case,  other  sr^eres  placed  at  e  =  ±  n/2  will  inter¬ 
fere  negligibly  with  the  original  sphere  since  the  probability  of  scat¬ 
tering  in  the  direction  e  =  ±  tt/2  is  almost  zero.  If  we  place  a  sphere 
along  8  =  0  we  will  certainly  be  faced  with  multiple  scattering  since 
the  probability  of  forward  scattering  is  large.  Thus  by  using  this 
simple  model  one  can  account  for  multiple  scattering  at  least  to  the 
first  order.  One  thing  we  have  not  discussed  is  the  phase  dependence 
on  the  position  of  the  scatterers  relative  to  each  other.  This  mutual 
dependence  leads  to  the  phenomenon  of  interference.  If  two  beams  of 
light  originate  from  the  same  source,  the  fluctuations  in  the  two  beams 
are  generally  correlated  and  the  beams  are  said  to  be  either  partially 
or  completely  coherent  [7],  If  two  beams  originate  from  two  independent 


sources,  they  are  said  to  be  mutually  incoherent  or  uncorrelated.  We 
can  therefore  conclude  that  for  scattering  by  two  spheres  placed  u't 
e  *  0  in  figure  1,  the  waves  scatte^d  by  the  two  spheres  will  be  at 
least  partially  coherent  since  the  incident  plane  wave  originated  from 
a  single  source.  We  can  also  consider  that  the  two  spheres  are  two 
different  sources  which  are  partially  correlated.  If  we  consider  a 
large  number  of  spheres  .placed  at  random  relative  to  each  other,  the 
scattered  waves  due  to  each  sphere  in  the  random  collection  will  be 
uncorrelated  and  fhe  scattering  will  thus  be  incoherent.  This  situa¬ 
tion  is  similar  to  the  problem  of  scattering  by  two  spheres'  that  are 
vibrating  randomly,  i.e.,  the  waves  scattered  from  them  will  be  uncor¬ 
related.  Thus,  if  we  assume  that  the  atmosphere,  for  example,  is  made 
up  of  many  randomly  placed  spheres*  we  can  determine  the  scattered 
intensity  at  any  point  incoherently,  i.e.,  without  regard  to  ohase. 

This  is  the  manner  in  which  we  will  approach  the  problem  of  scattering 
in  the  atmosphere. 

Depolarization  of  electromagnetic  waves  is  also  caused  by  aniso¬ 
tropic  media  where  one  or  more  of  the  characteristics  of  the  medium 
depends  on  direction.  In  anisotropic  media,  depolarization  takes  place 
because  of  the  reradiation  of  secondary  waves  with  characteristic 
polarizations  depending  on  direction.  Crystals  are  a  good  example  of 
anisotropic  medium  in  which  the  preferred  directic,-.  is  determined  by 
the  direction  of  the  magnetic  field.  Due  to  the  anisotropy  of  the 
medium  charges  can  move  only  in  a  specific  direction  and  are  not  free 
to  move  in  the  direction  of  an  incident  field  with  general  polarization. 
Only  if  the  incident  field  is  parallel  to  these  characteristic  direc¬ 
tions  will  the  polarization  of  the  reradiated  field  be  the  same  as 
the  incident  field. 
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Faraday  rotation  is  another  phenomenon  which  causes,  the  depol  ariza- 
tion  of  an  electromagnetic  wave.  Michael  Faraday  in  1845  and  Fresnel 
in  1837  found  that  when  light  propagates  through  a  medium  in  the  direc¬ 
tion  of  magnetic  field  the  plane  of  polarization  rotates  with  the  angle 
of  rotation  proportional  to  the  path  length  between  two  points.  The 
exact  derivation  of  this  phenomenon  can  bei  accompli  shed  only  through 
the  use  of  quantum  mechanics.  This,  form  of  depolarization  occurs  in 
the  earth's  ionosphere  [8]  and  can  cause  problems  for  radar  tracking 
sys  terns . 

In  this  study  we  wish  to  investigate  the  depolarization  of  electro¬ 
magnetic  waves  in  the  atmosphere.  We  will  consider  wavelengths  small; 
compared  to  interparricle  distances.  We  will  begin,  in  section  II  by 
discussing  the  radiative  transfer  of  energy  in  the  atmosphere  and 
derive  the  appropriate  integral  equations.  Section  II  follows  the  work 
of  Chandrasekhar  [9]  and  Kourganoff  [10]  closely.  The  class  of  integral 
equations  which  describe  the  problem  of  radiative  transfer  are  the 
nonhomogenecus  Fredholm  equations 

b 

viz)  =  F(0  *  *  J  Kfejx)  iU)  dx  (2) 

a 

which  has  been  solved  by  Chandrasekhar  [9]  for  the  isotropic  scattering 
function.  Evans  et  al.  [11]  have  produced  solutions  for  this  equation 
for  anisotropic  scattering  functions.  Both  Chandrasekhar  and  Chu  have 
used  Gauss  quadrature  in  solving  the  resulting  integral  equations.  The 
scattering  function  is  expanded  in  the  form 

K(n,n")  =  f(y) 

N 
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where 


m  =  cos  e 

e  =  angle  of  scatter 
=  angle  between  n  and  $r 
P.(y)  are  the  Legendre  polynomials 
and.  a |  are  the  coefficients  of  the  expansion.  Equation  (3)  has  the 
same  form  as  the  exact  expression  for  spheres. 

Solutions  of  equation  (2)  by  Gauss  quadrature  are  very  complicated 
and:  difficult  to  obtain.  In  this  study  we  obtain  numerical  results  to 
the  problem  of  depolarization  of  electromagnetic  waves  by  the  atmos- 
here  through  the  use  of  Monte  Carlo  techniques  [12„  13].  We  assume  that 
the  atmosphere  is  composed  of  spherical  scattering  particles  such  as 
water  droplets  in  clouds  or  fogs.  An  excellent  summary  of  various 
forms  of  atmospheric  conditions  is  contained  in  [1], 

in  section  ill  we  discuss  the  polarization  and  depolarization  of 
electromagnetic  waves  by  scattering.  We  also  present  the  Stokes 
parameters  in  describing  the  polarization  of  an  electromagnetic  wave. 
Section  III  follows.  Van  De  Hulst  [14]  closely  in  obtaining  an  analyti¬ 
cal  description  of  the  polarization  parameters.  Two  basic  forms  of 
scattering  by  particles  will  be  used  in  this  study.  One  is  Rayleigh 
scattering. where  the  scatterer  is  small  compared  to  a  wavelength  and 
the  other  is  Mie  scattering  where  we  consider  scattering  by  arbitrary 
spheres. 

Also  in  section  III  we  present  a  general  method,  as  was  obtained 
by  Deirmendjian  [15],  for  describing  the  size  distribution  of  particles 
in  the  atmosphere.  The  exact  solution  of  scattering  by  spheres  origin¬ 
ally  derived  by  Gustav  Mie  [2]  as  a  function  of  the  size  parameter 


will  be  used.  In  equation  (4)  r  is  the  sphere's  radius  and  X  is  the 
wavelength  of  the  incident  wave.  The  Mie  scattering  function  is  aver¬ 
aged  ov^er  the  size  distribution,  which  results  in  a  volume  scattering 
function.  If  f(e,r/x)  is  the  scattering  function- as  a  function  of  e, 
the  scattering  angle  and  the  ratio  r/x  is  as  defined  above;  and  n(r) 
is  the  si:ze  distribution,  then 

t 

F(e)  =  I  f(e *r/x)  n(r)  dr  (5) 

■  ri 

is  the  volume  scattering  function  for  particles  of  radius  between  rx 

and  r  . 

2 

In  section  IV  we  describe  the  computer  program  used  .for  simulating 
the  depolarization,  of  polarized  electromagnetic  waves  in  the  atmosphere. 
The  computer  program  is  based  on  previous  work  by  Kattawar  and  PI  ass  [;]6<j 
and  by  Collins  and  Wells  [17].  The  program  was  used  on  the  Air  Force 
Weapons  Laboratory  CDC-6600  computers. 

In  section  V  we  present  the  results  of  some  sample  problems.  We 
will  present  the  microscopic  scattering  functions  for  Mie  Scattering,  by 
"Cumulus  Clouds"  and  "Haze  M"  size  distributions.  By  microscopic  we 
mean  that  several  representative  size  particles  from  each  type  distri¬ 
bution  are  obtained.  We  will  also  obtain  and  present  the  volume  scat¬ 
tering  functions  for  each  size  distribution.  The  microscopic  and  volume 
scattering  functions  will  be  presented  for  several  different  wavelengths 
varying  from  0,5  micron  to  10.0  microns.  The  volume  scattering  functions 
are  th3n  used  to  analyse  different  atmospheric  problems  containing  either 
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the  "Cumulus  Cloud"  or  the  "Haze  M"  size  distribution  or  combinations 
of  both.  The  computer-program  has  been  written  to  accomodate  linear 
or  circular  polarization.  In  general,  the  distances  being  considered 
here  are  small,  i.e. ,  in  the  order  of  10  km  to  20  km  maxi  Ji.  Scat¬ 
tering  from  the  ground  is  included  and' we  will  therefore  include  the 
ground  albedo  in  each  problem.  The  albedo  is  defined  as  the  ratio  of 
scattered  to  incident  energy. 

The  result  of  all  these  problems  is  the  determination  of  the 
degree  of  depolarization  of  the  initial  source  electromagnetic  wave. 
Depolarization  for  the  vari qus  geometries  being  considered  is  pre¬ 
sented  for  different  wavelengths.  The  polarization  factor  defined  by 


D  = 


k  -  h 

‘r  +  I* 


(6) 


is  given  as  a  function  of  polar  angle;  where  Ip  is  the  perpendicular 
intensity  component  and  1^  is  the  parallel  intensity  component.  In 
equation  (6)  both  I  and  I  are  referred  to  a  vertical  plane  which 

I  X» 

contains  the  propagation  vector  (figure  2). 

Z 


Figure  2.  Vertical  Polarization  Reference  Plane 
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SECTION  II 


ELECTROMAGNETIC  SCATTERING  IN  THE  ATMOSPHERE 
1 .  INTRODUCf ION 

For  the  past  seventy  years,  astrophysicists  have  been  concerned 
with  the  transfer  of  radiation  in  the  atmospheres  of  planets  and 
stars.  Physicists  have  been  concerned' With  similar  problems  in  the 
diffusion  of  neutrons  in  various  materials.  We  should  also  mention 
that  radar  has  been  used- extensively  in  the  study  of  certain  atmos¬ 
pheric  properties. 

Recently,  with  the  invention  of  the  laser,  problems  of  light 
transfer  in  the  atmosphere  have  become  of  interest  to  a  greater 
number  of  scientists  and  engineers.  In  this  section  we  are  concerned 
with  the  formulation  of  the  problem  of  the  transfer  of  polarized 
light  in  the  atmosphere.  We  are  particularly  interested  in  the 
depolarization  effects  of  scattering  on  the  propagating  light  beam. 
The  notation  used  in  this  study  1s:;.  for  the  most  part,  that  used  by 
astrophysicists  and  follows  along  the  same  lines  as  in  the  books  by 
Chandrasekhar  and  Kourganoff  [9,  10], 

Basically,  the  problem  of  radiative  transfer  must  deal  with  a 
spherical  geometry  (if  we  can  assume  a  spherical  earth);  but  the 
problem  is  simplified  somewhat  if  the  atmosphere  is  plane  parallel. 

In  this  study  we  are  concerned  with  problems  with  short  range  trans¬ 
mission  of  highly  directive,  monochromatic  light,  i.e.,  the  distances 
involved  are  short  compared  to  the  curvature  of  the  earth.  It  is  also 
assumed  that,  at  the  frequencies  considered,  absorption  is  negligible. 
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2.  DERIVATION  OF  BASIC  QUANTITIES 

In  this  subsection  v:e  are  concerned  with  the  derivation  of  the 
essential  quantities,  necessary  for  the  solution  of  the  problem  as 
defined.  Since  all  quantities  defined  in  this  study  are  for  mono¬ 
chromatic  light  no  special  notation  is  used  to  denote  frequency 
dependence,  although  the  reader  should  keep  in  mind  that  all  quan¬ 
tities  are  highly  dependent  functions  of  frequency. 

Let  dE  be  the  energy  transported  across  an  element  of  area  dA, 
in  time  dt,  and  in  a  direction  confined  to  an  element  of  solid 
angle  dw  and  making  an  angle  0  with  the  normal  n  to  dA  (figure  3). 


dto 


n 

Figure  3.  Energy  Across  dA  in  the  Direction  dto. 

Then,  dE  can  be  expressed  in  terms  of  the  specific  intensity  or 
more  L  iefly  the  intensity  by 

dE  =  I  cos  e  dA  du>  dt  (7) 

If  the  intensity  I  is  independent  of  the  direction  at  a  point  the  radia 
tion  field  is  said  to  be  isotropic.  The  radiation  field  is  said  to  be 
homogeneous  if  the  intensity  is  the  same  at  all  points  and  in  all 
directions. 
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Radiation: which  is  propagating  through  a  medium  may  be  partially' 
scattered  or  absorbed  by  the  material  in  the  medium  which  in  essence 
weakens  or  attenuates  the  incident  intensity  I.  The  incident  inten¬ 
sity  I  becomes  I  +  dl  after  traversing  a  thickness  ds  .(figure  4)  in 
the  direction  of  propagation,  therefore, 

dl  =  -kp  I  ds  (8) 

where  k  is  the  mass  absorption  or  scattering  coefficient  and  p  is  the 
density  of  the  material.  In  what  follows  k  will  be  assumed  to  be  a 
scattering  coefficient:  exclusively. 


< 

•o 

- ►  dl 


Figure  4,  Intensity  Through;  Element  ds. 

Consider  a  pencil  beam  of  radiation  incident  on  an  element  of  mass 
with  cross-section  dA  and  thickness  ds.  The  scattered  intensity  in 
a  direction  dw  is  given  by 

dEg  =  («p  Ids)  dA  dw  dt  (9) 

The  mass  of  the  element  is 

dm  =  p  dA  ds  (10) 


%  * 

r4 

\ 

\ 


) 


Now  we  can  write 


I  dm  dw  dt 


4tt 


OD 
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Obviously,  in  the  most  general  case  the  scattered  energy  will  be 
dependent  on  direction;,  therefore,  a  scattering  function  F)(cos  o)  is 
introduced  to  account  for  this  dependence.  The  energy  scattered 
into  aa  element  of  solid  angle  dor*  making  an  angle  o  with  the  direc¬ 
tion  of  incidence  is  then 

dE$  =  k  I  P(cos  0)'  ^  dm  da.  dt  (12) 

and.  the  total  scs/ttered  energy  is 

E$  =  k  T  dm  dt  doi  P(cos  o)  doj^  (13) 

where  the  integration  is  carried  out  over  all  solid  angles.  In 
general,  for  perfect  scattering  (no  absorption) 

^  J P(cos  0)  dio'  =1  (14) 

in  general,,  where  both  scattering  and  absorption  occur 

V^gos  0)  *'<1  (15) 

where  is  the  fraction  cf  intensity  lost  due  to  scattering  and 
1  -  Uq  is  the  fraction  absorbed;  is  referred  to  as  the  albedo  for 
single  scattering.  If 

constant  =  P(cos  0)  06) 

the  scattered  intensity  is  isotropic. 

An  important  quantity  is  the  emission  coefficient  j.  The  emission 
coefficient  indicate0  the  amount  of  energy  emitted  by  an  element  of 
mass  dm  in  a  direction  of  solid  angle  do..  The  emitted  energy  includes 
scattered  energy  both  from  the  element  being  considered  and  any  other 
scattered  energy  originating  at  some  other  source  contained  in  the 
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solid  angle  d&>  being  considered.  Thus  the  energy  emitted  by  mass  dm 
in  the  direction  dw,  is  given  by 


Ee  =  j  dm  du  dt  (17.) 

the  scattering;  of  energy  incident  from  a  direction  (0",$")  'into  a 
di rection  -( e  ,<)>)  resul  ts  i n 

E-  =  ^  k  dm  dco  P(e,<j>;e",<j>'‘)  1(6-, <j> ")  sin  e"  de "  d4> "  dt  (18) 
Thus,  comparing  (17)  and  (18) 


2tt  ’tt 


j  =  ^  I  I  P(M*»e'><r)  I(e%<T)  sin  e"  de"  c 


0  0 


for  a  scattering  medium.  In  this  study  we  consider  scattering  only. 

The  ratio,  of  the  emission  to  the  absorption  coefficient  is  called 
the  source  function  and  is  written  as 


or  from  equation  '(19) 


2'd  IT 


0(e,$)  = 


‘fc/i 


P(e,^:e",<j>")  1(0", <j>")  sin  e'd  e"d  <t>' 


0  0 


3.  THE  EQUATION  OF  TRANSFER 

Consider  a  cylindericel  element  (figure  5)  of  cross-section  dA  and 
length  ds.  The  amount  of  energy  absorbed  by  the  element  is 


E,  =  ( tcp  I  d  s),  d  A  doj  dt 

a 


the  amount  emitted  is 
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dA 

Figure  5.  Cylindrical  Element  for  Derivation  of 
the  Equation  of  Transfer 

Eg  =  j  p  d  s  d  A  dwdt  (23) 

The  difference  between  the  absorbed  and  the  emitted  energy  results  in 
the  net  energy  crossing  the  two  faces  of  the  element  in  the  normal 
direction,,  and 


dl  dA  dwdt  =  E  -  E 
e  a 

or  using  (22)  and  (23) 

=  -kp  I  +  jp  (24) 

Using  equation  (20)  and  rearranging  we  obtain  the  equation  of  transfer 


dl 

<pds 


I  -  J 


(25) 


The  equation  of  transfer  (25)  can  be  solved  easily  through  the  use 
of  the  method  of  variation  of  parameters.  We  can  rewrite  equation 
(25)  as 


dl  ,  t  _ 

3?  +  pI  ~ 


<P 


Then 


I h ( s )  *  exP 


Kpd  s 


(26) 


(27) 
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where  0  is  a  constant  of  integration.  Combining  (27)  and  (33)  accord¬ 
ing  to  (28),  we  have 

s 

I  =  I(s)  =  C  e"7^s,°^  +  f  tcp  J(s')  e‘T^s,s^  ds“  (34) 

where 


t(s,sO 


(35) 


is  the  optical  thickness  of  the  material  between  s  and  s'.  At  s  =  0 
we  find  that 

1(0)  =  C 

therefore 


s 

I(s)'=  1(0)  e"x(s,°^  +  f  kp  j(s')  e"T^s,s^  ds'  (36) 

is  the  solution  to  the  equation  of  transfer  (25). 

For  the  case  of  scattering  in  a  plane-parallel  atmosphere  we 
measure  linear  distances  normal  to  the  plane  of  stratification,  then 
the  equation  of  transfer  becomes  (figure  6) 


t  =  0 


Figure  6.  Plane-Parallel  Atmosphere 
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-  cos  e  =  I(z,0  <f>)  -  J(z,0,$) 

were  the  angle  0  is  measured  from  the  z  axis  and  <p  is  the  azimuthal 


angle.  We  can  set 


dr 

3z  “  "Kp 


and  substitute  into  (37)  to  obtain 


y  d  ' dt  "  '  =  “  ^(T>Uj<i>) 


where  u  =  cos  0.  For  an  atmosphere  with  finite  optical  thickness  the 


solution  (36)  reduces  to  (figure  5) 


r( t ,+y ,<f» )  =  i(t,u,(J>)  e  (Vx)/lJ 


h 

■  I  J(t,y, 


*)  e" 


( t“T )/u  dt 


I  ( t  ,<f> )  -  I  (0 ,-u  ,(f>)  e  r//y 


•h 


where  1  >_  y  >  0  for  both  equations  (40)  and  (41).  Equation  (40) 
represents  the  intensity  in  the  positive  z  or  outward  direction, 
while  equation  (41)  represents  the  intensity  in  the  negative  z  or 
downward  direction.  The  intensities  emerging  from  the  atmosphere  at 
x  =  0  and  t  =  t  are  given  by 
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and 


I  (0  ,+y  ,<p)  =  I(t,p,<J>)  e  Tl^y 

+  /'  e-t/v  & 

11 


(42) 


^  (Ti  *‘J>)  =  1(0 »~m *<J>)  e  Ti//y 

f  J(tf-n,*)  e"(Ti"t)/,J  f-  (43) 

‘T) 

The  first  term  on  the  right  hand  side  of  equations  (40)  through  (43) 
may  be  interpreted  as  the  direct  intensity  attenuated  exponentially; 
the  second  term  is  the  sum  of  the  intensities  due  to  all  orders  of 
scattering  within  the  atmosphere. 

The  source  function  (equation  (21))  may  now  be  written  as 

2  71  +1 

j  (t,m)  ■  k/  /  p(»j><}>;p'><r) 

Jo  -1 

T(T,y%<r)  dp'df  (44) 


Writing 


Is(t,+m)  =  j  J  (t,y,*)  ™ 


(45) 


and 

T 

Is(t.-k,*)  *  f  0(t,-p ,♦)  e-l1-4)/"  ££  (46) 

■'o 

for  the  scatfar  contributions  in  equations  (40)  and  (41)  and  substi¬ 
tuting  eouation  (44)  into  (45)  and  (46),  we  obtain 
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P(n,<j»;y%<r) 


Tj  2lT 

h  =  hf  £ 

I(t, v',*')  d  u' d  ♦'  e^t'T'/“  ;p  (47) 
^and 

t  2tt  +1 

Ig (x  ,-M  ,«J>)  =  f  J  P{-W>^»u',<|>0 

I(t,M%r)  d  g'.d  r  er^’t)/v  ^  (48) 

If  we  now  compare  equations  (47)  and  (48)  with  equations  (40)  and  (41) 
we  can  see  that  we  now  have  a  set  of  integral  equations  to  solve  in 
order  to  acquire  the  intensity  at  any  point.  We  can  also  write 

Id(t>+M)  =  I  (tj  ,y  ,4>)  e-(TfTr  (49) 

and 

I[j(t,-u,<}»)  =  I(0,-u,<J>)  e  (50) 

as  the  direct  contributions  of  the  incident  intensity.  Therefore, 
using  equations  (47)  through  (50),  we  can  now  write  equations  (40)  and 
(41)  as 


I  (  T  ,+M  ,4>  )  =  Ip^Tjj+U,^  +  Is(t,+W,<J>) 


(51) 


and 


I(t,“u,4>)  =  I|j(t*-m>4>)  +  Is  (x  ,4>) 


(52) 


i.e.,  the  intensities  propagating  in  positive  and  negative  directions 
can  be  considered  as  being  made  up  of  a  direct  contribution  attenuated 
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exponentially  and  a  contribution  consisting  of  the  intensity  from  all 
orders  of  scattering. 

4.  RADIATIVE  TRANSFER  .FOR  POLARIZED  LIGHT 

Polarized  light  can  be  described  in  terms  of  the  Stoker  Dara.neters 
(section  III).  There  are  four  Stokes  .parameters  which  describe  the 
total  intensity,  the  polarization  ellipse,  and  the  orientation  of  the 
polarization  ellipse;  The  four  Stokes  parameters  are  written  as 

1=  (I„Q,U,V)  (53) 

where  the  boldface  letters  will  indicate  a  matrix  representation. 

Procedi ng  with  matrix  notation,  we  can  now  represent  the  scattering 
function  by  a  scattering  transformation  matrix.  The  transformation 
matrix  will  be  .a  4  X  4  matrix  and  is.  of  the  form  (see  Section  III)  [18] 


/31  bl  b3  b5\ 


C3  C4  a3  k2 


\C5  C6  C2  a4/ 

Furthermore,  if  we  write  the  polarization  parameters  as 

x  *  ( Wu,v)  (56> 

then  the  matrix  transformation  will  be  denoted  by  F*.  We  will  find, 
in  section  III,  that,  for  a  spherical  scattering;  particle,  the  scat¬ 
tering  matrix  will  be  of  the  form 


22 


F  = 


/a  {  0 


a2 


\0  0 


0 

o- 


a3 


0 

0 


b2 

a4, 


(56) 


where  b;  =  -c;  and  a"  =  a,". 

2  2  3  4 


Using  matrix  notation  we  can  now  write  equations  (51)  and  (52)  as 


and 


where 


I  (  T  ,+P  »<}> )  = 

Iq(t  ,+y  ,<J>)  +  I$(x  ,+p  ,4>) 

(57) 

I  ( T ,-u  ,4>)  = 

ID(  t  ,-p  ,<f>)  +  Is(t,-u,*) 

(58) 

Iq(t,+U,<J>) 

■  l(r,  e-(TrT)/lJ 

(59) 

In(T,-u,<i>) 

=  I(0,-u,<fr)  e"T/lJ 

(60) 

1 


T,  2ir  +1 

/// 


Is(t,+u,<j>)  =  I  J  J  F  (iMlu',^) 

0  -1' 


Kt.u'.r)  d  u' d  *'  e-^-T),u  J1 


(61) 


T  2TT  +1 

Ig ( t »"U )  =  ^  J"  J*  F  (“Wsi’jW  ) 

0  0  -1 

I(t,p',r)du'df  e'^f  (62) 

Equations  (57)  through  (62)  describe  the  intensity  of  radiation  at  any 
point  in  a  semi -infinite,  plane  parallel  scattering  atmosphere. 

The  solutions  (57)  and  (58)  are  nonhomogeneous  Fredholm  integral 
equations.  The  difficulty  in  solving  this  type  of  integral  equation 
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is  that  the  function  being  solved  for  also  appears  under  the  integral 
equation.  Generalized  physically  realizable  kernals  may  be  solved 
numerically  through  the  judicious  use  of  series  expansions. 

5.  MODELING  THE  TRANSFER  PROBLEM 

A  study  of  equations  (57)  and  (58)  will  help  in  envisioning  the 
physical  events  which  take  place  during  the  propagation  of  energy 
through  a  scattering  atmosphere.  First,  as  already  indicated,  the 
intensity  at  any  point  can  be  thought  of  consisting  of  the  direct 
intensity  attenuated  exponentially  and  of  the  scattered  intensity  from 
all  scatterers  surrounding  the  observation  point.  The  polarization 
properties  of  the  direct  beam  do  not  change,  while  the  scattered 
intensity  will  in  general  have  its  polarization  properties  modified  by 
the  scattering  function  matrix  of  each  individual  scatterer. 

The  intensity  of  electromagnetic  radiation  is  defined  [19]  as  the 
energy  per  unit  time,  passing  through  a  unit  area  at  right  angles  to 
the  direction  of  propagation.  In  terms  of  a  beam  of  monochromatic 
photons  the  intensity  is  given  by 

I  =  Ep  (Photon  Flux)  (63) 

where  E  =  hv,  h  is  Planck's  constant  and  v  is  the  photon  frequency. 

r 

Equation  (63)  suggests  that  a  beam  of  monochromatic  photons  is  the 
proper  representation  for  directed  electromagnetic  radiation.  Suppose 
we  approximate  equations  (61)  and  (62)  by  the  triple  summation  of  a 
large  number  of  terms 

n  m  l 

W  e’(T 

i=l  j=l  k=l 
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this  triple  summation  is  equal  to  (61)  in  the  limit  as  A»T,A<J>',  and 
At  go  to  zero  and  n,  m,  i  go  to  infinity.  Furthermore,  the  beam;  may 
be  thought  of  as  being  composed  of  many  photons  according  to  equation 
(63).  We  may  therefore  simulate  a  beam, -of  radiation  by  individual 
photons  or  individual  photons  individually  propagating  and  scattered 
randomly  about  throughout  the  medium.  Since  we  are  only  considering 
noncoherent  scattering,  the  results  may  be  added  arithmetically  at 
some  point  of  observation.  It  will  further  be  assumed  that  near 
field  effects  are  negligible. 

The  essence  of  the  model  is  this:  Consider  a  photon  or  group  of 
photons  with  a  defined  set  of  Stokes  polarization  parameters,  launched 
from  an  isotropic  source;  these  photons  will  be  attenuated  exponen¬ 
tially  as  they  propagate  through  the  medium,  and  upon -.colliding  with 
a  particle  the  photons'  polarization  parameters  will  be  modified 
according  to  the  scattering  matrix  and  scattered  in -a  new  direction, 
the  history  of  this  group  of  photons  may  be  followed  to  all  orders  of 
scattering.  In  this  .way  many  histories  may  be  followed  and  estimates 
of  the  intensities  at  a  detector  (observation  point)  may  be  made. 

This  technique  has  been  used  by  other  investigators  [16,  17,  20]  with 
excellent  agreement  with  experimental  results  for  natural,  i.e., 
unpolarized  electromagretic  radiation.  In  this  study  we  use  this 
approach  for  the  case  of  a  polarized  light  source  and  determine  the 
effect  of  scattering  on  the  initial  polarization  state  of  the  beam 
of  light. 

In  the  next  section  we  discuss  the  Stokes  parameters  and  deter¬ 
mine  the  effect  of  scattering  on  these  parameters. 
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SECTION  III 


POLARIZATION  AND  DEPOLARIZATION  OF  EM  WAVES 
1 .  INTRODUCTION 

Electromagnetic  waves  have  several  properties  which  can  be  used 
to  convey  information  from  one  point  to  another.  The  most  commonly 
used  properties,  are  the  frequency,  phase,  and  amplitude.  The 
electromagnetic  wave  propagating  in  a  medium  whether  it  be  a  wire, 
a  waveguide,  or  the  atmosphere  can  be  used  to  carry  information  by 
modulating  the  frequency,  the  phase,  or  the  amplitude  according  to 
some  prescribed  code.  The  information  can  be  recovered  at  the 
receiving  end  through  demodulation. 

Another  property  of  electromagnetic  waves  which  has  not  been 

used  to  a  great  extent  is  the  polarization  of  the  wave  propagating 

through  a  medium.  It  is  known  that  the  electric  and  magnetic  field 

vectors  are  transverse  to  the  direction  of  propagation  in  an  isotropic, 

homogeneous  medium.  The  orientation  of  the  electric  field  vector 
*  > 

defines  the  polarization  of  the  electromagnetic  wave.  By  convention, 
if  the  direction  of  the  electric  vector  is  constant  in  time  the  wave 
is  said  to  be  linearly  polarized,  if  electric  vector  rotates  then  the 
wave  is  said  to  be  elliptically  polarized. 

The  polarization  of  an  electromagnetic  wave  propagating  through 
a  continuum  remains  unchanged  if  the  medium  is  homogeneous  and  iso¬ 
tropic.  On  the  other  hand,  if  scattering  occurs  the  polarization 
characteristics  will,  in  general,  change.  Therefore,  if  the  polari¬ 
zation  characteristics  of  a  wave  prior  to  and  after  scattering  are 
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known,  it  is  possible,  at  least  in  principle,  to  categorize  the 
particular  type  of  scatterer  causing  the  depolarization.  In  this 
study,  we  are  primarily  interested  in  the  de do! aeration  of  a  com¬ 
pletely  polarized  wave  by  a  scattering  atmosphere.  By  depolariza¬ 
tion  we  mean  the  degree  to  which  an  arbitrarily  polarized  transmitted 
wave  is  transformed  to  some  other  polarization  state. 

2.  POLARIZATION  AND  THE  STOKES  PARAMETERS 

The  most  general  type  of  polarized  wave  can  be  written  as 

E  =  \  +  rEr  (65) 

where  < 

•f 

E  is  the  total  electric  field  vector 
E£  is  the  parallel  component 
Er  is  the  perpendicular  component 
z  and  r  are  unit  vectors 

Jn  writing  equation  (65)  we  used  the  same  notation  used  by 
Chandrasekhar  [9]  and  Van  De  Hulst  [14],  where  the  subscripts  z 
and  r  refer  to  the  last  letters  in  the  words  parallel  and  per¬ 
pendicular,  respectively.  The  parallel  and  perpendicular 
components  are  referenced  to  some  arbitrary  plane,  usually  the 
plane  of  scattering.  In  general,  we  can  write  each  component  as 

Et  =  e'ci  e‘1kz  *  iwt  (66) 

and 

Er  =  ar  e~Uz  e~1kz  +  ia)t  (67) 

where 
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a0  and  a,,  are  the  positive  real  amplitudes 
£X  and  e2.  are  the  phase  angles 
k  =  2u/x  is  the  propagation  constant 
a  =  2-irf  is  the  monochromatic  angular  frequency 


Equations  (66)  and  (67)  represent  a  plane  wave  propagating  in 
the  positive  z  direction.  These  equations  also  represent  a  general 
elliptically  polarized  wave.  In  this  study  we  are  using  the  Stokes 


parameters  to  characterize  the  polarization.  The  Stokes  parameters 

(I,  Q,  U,  V)  are  defined  as 

1  ■  E*  **r  E* 

(68) 

0  ■  E*  EI  '  Er  K 

(69) 

U  =  E*  ^  +  EA  Er 

(70) 

v  ■  f(Et  k  -  $  Er) 

(71) 

from  (66)  and  (67)  we  have 

E*  EI  ■  n* ' 

(78) 

E  E*  =>  a  2 
r  r  r 

(73) 

c  c*  -  (ei  "  ^2) 

Er  “  aa  ar  e 

(74) 

r*r  i(*l  “  *2) 

E*  Er  =  h  ar  e  v 

(75) 

applying  equations  (72)  through  (75)  to  equation:*  (68)  through  (71) 
we  obtain 
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(76) 


where 


I  =  a.2  +  a  2 
z  r 


0  -  V  -  ar2 


U  =  2a.  a^  cos 6 
£  r 


V  =  2a,  a  sins 


(77) 

(78) 

(79) 


rrom  equations  (76)  through  (79)  we  see  that  for  parallel  polariza¬ 
tion  (ar  =  0) 


1  =  V 

Q  -  a/ 

U  =  0 

V  =  0 

for  perpendicular  polarization  (a^  =  0) 


U  =  0 
V  =  0 


.(80) 

(81) 

(82) 

(83) 

(84) 

(85) 

(86) 
(87) 
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For  circular  polarization,  we  must  define  a  convention  for  describ¬ 
ing  the  sense  of  rotation.  For  this  we  use  the  standard  convention 
which  defines  right  handed-  polarization  with  a  wave  that  is  viewed 
traveling  away  from  the  observer  and  whose  electric  vector  rotates 
in  the  same  direction  as  a  right  handed  screw.  Similarly,  left 
handed  polarization  is  defined  with  a  wave  viewed  traveling  away 
from  the  observer  and  whose  electric  vector  rotates  in  the  opposite 
direction  to  that  of  a  right  handed  screw.  Figures  7a  and  7b >show 
typical  right  hand  and  left  hand  polarization,  respectively. 

Now,  let  us  write 


E£  =  a-  cos(wt  -  6)  (88) 

Er  =  ar  cos  wt  (89) 

If  we  let  6  =  +  ti/2 ,  we  obtain 

=  a£  cos  (wt  -  u/2)  =  a^  sin  wt  (90) 

But,  since 

E  =  z  E.  +  r  -Ew  *  z  a.  sin  wt  +  r  a„  cos  wt  (91) 

z  r  z  *  r 

this  is  a  general  right  hand  elliptically  polarized  wave  (figure  8). 
From  figure  8  we  can  also  see  that  if  af  *  a^  right  nand  circular 
(RHC)  polarization  is  obtained.  If  6  =  -  tt/2,  equation  (88)  becomes 

E^  =  a£  cos  (wt  +  tt/2)  =  -  a£  sin  wt  (92) 

and 

E  =  -  s,  a&  sin  wt  +  r  ar  cos  wt  (93) 
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z 

direction  of 
propagation 


X 


direction  of 
propagation 


Typical  Right  and  Left  Hand  Polarization 


a 


Figure  8.  Right  Hand  Elliptical  Polarization 

Equation  (93)  is  the  general  representation  of  a  left  hand  ellipti 
cally  polarized  wave.  Figure  9  shows  equation  (93)  for  a  =  a  = 


Figure  9.  Circular  Polarization  for  ar  =  a£  =  aQ 

Therefore,  letting  <$  =  +  n/ 2  and  a&  =  ar  =  aQ,  the  Stokes  parameters 
for  right  hand  circular  polarization  are 


Q  =  O' 

(95) 

U  =  0 

(96) 

V  =  2a§ 

(.97) 

Similarly,  if  6  = 

-  tt/2  and  ar  =  a&  =  ao*  the  Stokes  parameters 

for 

left  hand  circular 

•  polarization  (LHC)  are 

I  =  2a§ 

(98) 

Q  =  0 

(99) 

U  =  0 

(100) 

V  •  -  2a§ 

(101) 

If  we  let 

I,  ■  at2 

and 

*r  *  ar  ' 

where  I.  and  I  are  the  intensities  for  the  parallel  and  perpendic- 
x»  r 

ular  intensities. 

respectively,  we  can  rewrite  equations  (76) 

through  (79)  as 

- 

-  -  -- 

'  I  =  I*  +  Ir 

(102) 

- 

Q  *  h  -  ir 

(103) 

U  =  2^ll  If  cos<5 

(104) 

V  =  2^  Ir  sin6 

(105) 
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A  more  intuitive  representation  of  the  polarization  would  be  to  use 
the  set  (I^,,  Ir,  U,  V)  rather  than  (I,  0,  U,  V).  This  new  set  of 
parameters  is  more  pleasing  in  that  the  relative  si/ze  of  parallel 
and  perpendicular  components  can  be  quickly  determined-. 

Up  to  this  point  we  have  discussed  the  Stokes  polarization 
parameters  for  the  ease  where  the  major  and  minor  axis  of  the 
polarization  ellipse  are  along  the  z  and  r  axes.  In  general,,  this 
may  not  be  true  although  a  simple  rotation  of  axis  would  realign 
the  axes  in  the  propei?  orientation.  Suppose  the  axis  of  the  .polari¬ 
zation  ellipse  were  rotated  by  an  angle  x  (figure  10) 
r 


r 


along  the  new  axes  p  and  q.  The  electric  vector  can  be  written  as 

E  =  a  p  cose  sin(wt-kz)  +  a  q  sine  cos(ut-kz)  (106) 

where 
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a  cose  is  the  amplitude  along  the  p  axis 
a  sine  is.  the  amplitude  along  the  q,  axis' 
p  and  q  are  unit  vectors 

If  we  let  e  ~  0  in  (106),  than 

E  =  a  p  sin(wt  -  kz)  (107) 

which  is  a  linearly  polarized  plane  wave.  Similarly  if  e  =  tt/4 
we  ;have  RHC  .polarization  and1  if  e  ==  -  tt/4  we  have  LHC  polarization. 

In  general,  the  ellipticity  is  given  by  tane,  if  it  is  0  we  have 
linear  polarization,  if  it  is  -1  we  have  LHC  polarization  and  if 
it  is  +1  we  have  RHC  polarization.  The  orientation  of  the  ellipse 
is  given  by  the  angle  x*  From  figure  10  we  have 


r  =  p  sin  x  +  q  cos  x 

(108) 

t,  =  p  cos  x  -  q  sin  x 

(IDS) 

Rewriting  equations  (66)  and  (67) 


Et  =  at  e"1*1  e-,kz  +  1ut 

((66) 

Er  .  ar  e'1'2  e'1kz  +  iBt 

(67) 

we  have  that 

E  «  Re  /e,  «  +  Er  r j 

(110) 

Substituting  (f 

and  (67)  into  (110)  results  in 

E  =  *«t 

cos  (wt  -  kz  -  e»)  +  r  ar  cos  (wt  -  kz  -  e2) 

(111) 
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Upon  equating  equations  (106)-  and  (111),  introducing  equations 
(108)  and  (109),  and  equating  Mke  terms  we  have 

a^  cos  (ut-kz-ej)  =  a  cos  8  cos  x  sin  (wt-kz) 

-  a  sin  e  sin  x  cos  (wt-kz)  (112) 

ar  cos  (iot-kz-e2)  =  a  cos  e  sin  x  sin  (wt-kz) 

+  a  sin  B-  cos  x  cos  (wt-kz)  (113) 


We  can  now  apply  the  relations 


A  cos  e  +  B  sin  e  =  \/a2+B2  cos  fc-tari"*  | ^ 


and 


A  sin  0  -  B  cos  e  =  y/kz+Q2  cos^e+tan"1  ^  j 
to  thev  right  hand  side  of  equations  (112)  and  (113)  as  follows 


(114) 

(115) 


afc  cos  (wt-kz-ej)  =  a  Vcos2  3  cos2  x  +  sin2  8  sin2  x 

cos  (uit-kz  +  tan"1  -c-?-s-  A  X,\.  (116) 

l  sin  8  sin  x/  '  ' 

ar  cos  (mt-kz-En)  =  a  Vcos2  6  sin2  x  +  sin2  8  cos2  x 

-  (-w*  -  *»-*  017) 

so  that 

a2  =  a2  (cos2  8  cos2  x  +  sin2  8  sin2  x)  (118) 

a2  =  a2  (cos2  8  sin2  x  +  sin2  8  cos2  x)  (119) 
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or  since 


tan  e:  =  -  cot  8  cot  x 

(120) 

tan  e£  =  cot  e  tan  x 

(121) 

Ej-e 2  =  <$  =  -tan_1(cot  8  cot  x)-tan_1(cot  8  tan  x)  (122) 
we  have  that 


tan 

tan'1 (cot  8  cot  X) 

+ 

tan. 

- 

stfaii- 1  (cot  8  tan  x) 

lrtan  tan"1 (cot  6  cot 

x) 

tan 

tan~1.(cot  8  tan  x) 

_  cot  8  cot  x  +  cot  8  tan  x 
cot2  -8  -  1 


1  tan2  X  +  1 
■  .  — —  ..,  ■■■■  «  . —  !■■■■ 

cot  2  8  2  tan  y 


-  tan  2  8 
”  sin  2  x 


(123) 


Using  equations  (118),  (119),  and  (123),  the  Stokes  parameters 
become 

(124) 

(125) 

To  obtain  cos  6  we  use  equation  (123)  and  figure  11. 


I  =  a2  +  af2  =  a' 


Q  ®  a|  -  ar2  =  a2  cos  2  8  cos  2  x 


U  =  2a.  a^  cos  6  =  2a2(cos2  8  cos2  x  +  sin2  6  sin2  x) 

X/  i 


1/2 


(cos2  6  sin2  x  +  sin2  8  cos2  x)cos  6 
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tan  26 


Figure.  11.  Triangle  for  Deriving  Cos  6 
Therefore  ,  we  have 


cos  6  = 


VT 


sin  2X _ 

an2  2b  +  sin2  2x 


and  using  straightforward  algebra 


U  =  a2  sin  2x  cos  2b 
Again  from  figure  il  we  obtain 


and 


sin  6 


tan  2B 

Vtan2  2B  +  sin2  2x 


(126) 


V  =  a2  sin  2b  (127) 

Equations  (124)  through  (:127)  lead  us  to  another  way  of  describing 
the  polarization  state  of  an  electromagnetic  wave,  namely  through 
the  use  of  the  Poincare  Sphere. 

For  convenience  we  now  rewrite  equations  (124)  through  (127) 


I  = 

a2 

(128) 

Q  = 

a2  cos  2b  cos  2x 

029) 

U  = 

a2  cos  2b  sin  2x 

(130) 

V  = 

a2  sin  2b 

031) 
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If  P  is  any  point  on  a  sphere  (see  figure  12)  it  can  be  described 
by  the  parameters  (I,  Q,  U,  V). 


Figure  12.  Representation  of  the  Polarization  Parameters 
with  the  Poincare  Sphere 

Comparing  equations  (124)  through  (127)  with  equations  (94)  through 
(97)  we  can  see  that  if  23  =  tr/2  we  obtain  right  hand  circular 
polarization;  this  is  the  north  pole  of  the  sphere.  Comparing 
equations  (128)  through  (131)  with  equations  (98)  through  (101)  and 
letting  2e  =  -  n/2  we  obtain  left  hand  circular  polarization;  this 
is  the  south  pole  of  the  sphere.  If  23  =  2x  =  0,  we  get  parallel 
polarization;  for  2s  =  0  and  2x  =  w,  we  get  perpendicular  polariza¬ 
tion.  In  general,  the  upper  hemisphere  yields  right  handed 
elliptical  polarization,  the  lower  hemisphere  yields  left  handed 
elliptical  polarization,  and  the  equator  yields  linear  polarization. 
The  total  intensity  is  given  by  the  magnitude  of  the  radius  vector. 
Shurcliff  [21]  discusses  the  graphical  mapping  of  one  polarization 
state  to  another  through  the  use  of  the  Poincare  Sphere. 
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3.  RESOLUTION  OF  A. GENERAL  ELLIPTICALLY  POLARIZED  WAVE  INTO  RHC 
AND  LHC  COMPONENTS 

According  to  the  previous  subsection,  an  jlliptically  polarized 
wave  may  always  be  resolved  into  a  perpendicular  polarization 
component  and  a  parallel  polarization  component  with  a  phase  dif¬ 
ference  between  them.  In  like  manner.,  an'  elliptically  polarized 
wave  can  be  resolved  into  an  RHC  component  and  an  LHC  component. 

This  latter  representation  yields  more  information  than  the  former 
because  the  phase  difference  is  automatically  included.  In  this 
subsection,  a  general  elliptically  polarized;--w.aye  will  be  resolved 
into  its  RHC  and  LHC  components. 

Consider  a  polarized  wave  of  the  form 

E  =  4  E*  +  r  Er  (65) 

as  before,  with 

E&  =  a£  cos  cut  (132) 

Er  =  ar  cos  (at  -  6)  (133) 


Each  component,  E  and  E  ,  may  be  individually  resolved  into  RHC  and 

X>  i 

LHC  components,  i.e., 


*  E£  “  EtRHC  +  EaLHC 


(134) 


with 


**■  +  . 
EARHC  ~  T"  ^  C0S  "  r  s^n 


(135) 


•tLHC 


.  a* 


(£  cos  wt  +  r  sin  u>t) 


(136) 
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While  for  the  perpendicular  component 


rEr  =  ErRHC  +  ErLHC 

(137) 

where 

ErRHC 

=  [ z  sin  (wt-6)  +  r  cos  (wt-6)] 

(138) 

ErLHC 

=  y~  [-t  sin  (wt-6)  +  r  cos  (wt-6)J 

(139), 

Therefore,  from  equations  (65),  (134),  and  (137) 

-> 

E  =  i 

(E*RHC  +  E£LHc)  +  (ErRHC  +  ErLHc) 

(140) 

or 

M 

^RHC  +  ErRHc)  +  (EaLHC  +  ErLHc) 

(141) 

letting 

r>*  4 

ERHC  =  E£RHC  *  ErLHC 

(142) 

and 

-y  -y  -y 

ELHC  =  EALHC  +  ErLHC 

(143) 

we  arrive  at 

E  =  ERHC  +  elhc 

(144) 

Substituting  equations  (135),  (136),  (138),  and  (139)  into  (144) 
and  writing  the  RHC  and  LHC  components  separately  we  have 
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•RHC  '  2 


0  ^ 

-  a  cos  wt  -  r  sin  wt 


+  if-  jt  sin  (wt  -  6}  +•  r  cos  (wt  r  6)J 


=  1/2  t  |a&  cos  wt  +  ar  sin  (wt  -  6) 


1/2  r  -a.  sin  wt  +  a^  cos  (wt  -  <$•) 

~  *  Jc 


(145) 


ao  r*  ->■  n 

ELHC  =  2“  I .  cos  +  r  sin  wtj 


+  2^"  [-1  sin  (wt  -  6)  +  r  cos  (wt  -  $)J 


=  1/2  t  ^  cos  wt  -  arsin  (wt  -  &)J 
+  1  2  r  ja£  sin  wt  +  afcos  (wt  -  6 )J 


For  (145)  and- (146)  we  can  write 


(146) 


sin  (wt  -  6)  =  sin  wt  cos  6  -  cos  wt  sin  6 


(147) 


COS  (wt  -  6)  =  COS  wt  cos  6  +  sin  wt  sin  6 
Substituting  (147)  and  (148)  into  (145)  and  (146)  we  obtain 


(148) 


Erhc  =  X  j^(a£  -  af  sin  6)  cos  wt  +  ar  cos  6  sin  wtj 

+  1/2  r  (az  -  af  sin  6)  sin  wt  +  ar  cos  6  cos  wtj  (149) 
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and 


ELHC  =  ^2  *  £(a£  +  af  sin  6)  cos  wt  -  ar  cos  <5  sin  wtj 

+  1/2  r  [(a£  +  ar,  sin  s)  sin  wt  +  ar  cos  6  cos  wtj  (150) 


Now,  equations  (149)  and  (450)  are  of  the  form 


A  cos  wt  +  B  sin  at  =  y^A2  +  B2  cos  ^wt  -  tan"1 


and  may\be  written  in  this  form  as  follows: 


-RHC  =  1/2  *  |AaV  "  ar  s1n  6)  +  (ar  cos  6) ' 


(a  cos  <$ 

wt  -  tan"1  — — - 

a£  -  ar  sin  6 


(151) 


+  1/2  r|/(a£  -  ar  sin  <$)2  +  (ar  cos  6): 


cos  f, wt  +  tan 


_j  a£  -  ar  sin  6 


a„  cos  6 
r 


(152) 


E,ur  *  1/2  a -f/Ca.  +  a  sin  s)2  +  (a„  cos  6)2 


'  '  x,  y  v  / 

(a  cos  6 

wt  +  tan-1  — - - 

a£  ar  sin  6 


+  V2  r  |/(a£  +  ar  sin  6)  +  (ap  cos  6)2 


(a  +  a  sin  6 

wt  -  tan"1  — - - - 

ar  cos  6 


(153) 


Furthermore  we  may  write  from  equation  (152) 
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^RHC 


aM  cos  6 
=  tan-1  — T— — - 


^rRHC  =  tan 


a^  -  ar  sin  6 


a„  -  a„  sin  6 


-l  _£ _ 1 


aM  cos  6 

r 


with 


(see  figure  13). 


^RHC  =  ff/2  “  * 


rRHC 


arcos6 


Figure  13.  Triangle  for  Deriving  <j>rRHC  and  <f^RHC 
Similarly,  from  equation  (153) 


**LHC  =  tan 


a„  cos  6 
■l  r 

a 5i,  '+  ar  sin  6 


♦rLHC  “  tan 


a,  +  a„  sin  s 
_  i  a  r _ 

ar  cos  6 


with 


(-154) 


(155) 


(156) 


(157) 


(158) 


*ILHC  =  "/2  '  * 


rLHC 


(159) 


(see  figure  14). 
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Figure  14.  Triangle  for  Deriving  <j>r^c  and  <^HC 

Using  equations  (154)  through  (159)  and  carrying  out  the  squaring 
of  the  terms  under  the  radicals  equations  (.152)  and  (153)  become 

Erhc  *  1/2  Va*  ■  2as  ar  s1n  s  +  *1 

t  cos  (»t  -  *tRHC)  +  r  cos  (»t  -  *1RHC  +  ir/2)J 


but 


?LHC  =  1/2  Va*  +  2at  ar  s1n  6  +  a? 

[*  cos  (“‘  +  +sLHc)  +  ?  cos  (“‘  +  +SLHC  -  "/2)] 
cos  («t>  ♦JtRHC  +  «/z)  ■  -  sin  -  *tRHc) 


and 


cos  (»t  +  otLHC  -  s/2)  ■  +  sin  (at  +  *tLHC) 


therefore, 


Er„c  =  1/2  ya|  +  2a*  0r  sin  6  +  a* 

[*  cos  (“'  -  *Mc)  '  ?  s1n  (st  ’  +*RHc)J 


(160) 
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(161) 


!lhc  -  1/2  +  2at'sr  stn  «  +  a* 

[‘  cos  (“‘  +  **LHC)  +  ?  s1n  (“*  +  *tLHc)] 


The  phase  difference  is 


*  =  **RHC  +  %HC 


a„  cos  6  h  cos  6 

tan-1 - +  tan"1  — - - « — 

a„ -aM  sin  6  a.  +  a  sin  6 

*  r  £  r 


which,  when  using  the  identity 


becomes 


tan"1  x  +  tan"1  y  = 


<s>  =  tan"1 


(162) 


Equations  (*160)  through  (162)  constitute  the  RHC  and  LHC  components 
of  a  general  ellvptically  polarized  wave.  Note  that  <t>  is  defined 
only  when  the  wave  is  not  circularly  polarized.  When  the  wave  is 
circularly  polarized,  either  the  RHC  or  LHC  component  will  be  zero 
and  $  will  not  be  defined.  That  this  is  true  car.  be  seen  by  letting 
a£  =  ar  and  6  =  +  n/2  simultaneously  as  is  required  for  circular 
polarization. 

4.  DEPOLARIZATION  AND  POLARIZATION  FACTORS 

In  this  study,  depolarization  will  be  defined  as  the  change  in 
polarization  state  which  includes  the  completely  depolarized  state  as 
natural  light.  The  definition  f  natural  or  unpolarized  light  is 
taken  from  Chandrasekhar  [9]: 
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Natural  Tight  is  equivalent  to  any  two  independent 
oppositely  polarized  streams  of , half  the  intensity-^ 

.ancLnc  two  independent  polarized  streams  can  to¬ 
gether  be  equivalent  to  natural  light  unless  they 
be  oppositely  polarized  and  of  equal  intensity. 

Experimentally,  if  natural  light  is  resolved  in  any  direction  in  the 

transverse  plane,,  the  measured  intensity  is  the  same.  The  necessary. 

and  sufficient  condition  for  light  to  be  natural  is  that 

q  =  u  =  V  =  0  (163) 

A  theorm  due  to  Stokes  [22]  states  that  the  most  general  mixture 
of  light  can  be  regarded  as  a  mixture  :of  an  elliptically  polarized 
wave  and  an  independent  wave  of  natural  light. 

From  equations  (68)  through  (71)  we  can  show  that 

12  =  q2  +  (J2  +  v2  (164) 

But,  for  the  most  general  representation  we  can  write 

I2>q2  +  u2  +  v2  (165) 

since  Stokes  theorem  says  we  can  write  the  Stokes  'parameters  as 

I  ■  I,  +  Ir  +  *„  (’66) 

Q  *  I,  -  Ir  (167) 

U  =  2  VTT’r  cos  6  (168) 

V  =  2  V’s.  'r  sin  4  069) 

Where  I  is  the  intensity  for  natural  light.  The  degree  of  polariza¬ 
tion  of  a  wave  is  defined  as  [9] 
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(170) 


n  _  \^Q2  +  U2  +.  V-2 
K  f 

and  from  equation  (165)  it  follows  that 

0<P<1  (171).- 

When  P  =  1  the  wave  is  said  to  be  completely  polarized)  when  P  =  0 
the  wave  is  said  to  be  completely  unpolarized  (natural  light).  In 
general  _  .  .  - 


0  <  P  <  1 


(172) 


and  the  wave  is  partially  polarized. 

In  this  study,  we  will  be  dealing  with  the  transmission  of  com¬ 
pletely  polarized  (P  =  1)  or  completely  unpolarized  light  through 
model  atmospheres.  We  will  be  interested  in  the  affect  of  scattering 
on  the  polarization  properties  of  the  transmitted  wave.  Comparisons 
with  natural  or  unpolarized  light  will  be  made.  It  will  further  be 
assumed  that  due  t  the  randomness  of  the  particle  distribution  in 
the  atmosphere,  the  intensities  of  the  individual  scattering  events 
can  be  summed  at  the  detector,  the  transmitted  energy  will  be 
monochromatic,  and  the  individual  scattering  particles  are  located 
far  enough  from  each  other  such  that  the  far  zone  approxi'  .tion  may 
be  used. 

Since  we  are  interested  in  the  degree  of  depolarization,  i.e., 
the  amount  of  cross-polarized  component  produced  by  scattering  of 
a  polarized  wave  in  the  atmosphere,  therefore,  in  an  effort  to 
provide  qualitative  results  and  for  comparison  of  various  cases,  we 
will  introduce  the  polarization  factors.  For  linearly  polarized  and 
unpolarized  transmitted  waves,  we  will  use 
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(173) 


lr  i  lZ 
D,  =  — - ~ 

V  +  I£ 


while  for  circularly  polarized  transmitted  waves  we  will  use 


D 


C 


=  *RHC  "  !LHC 
*RHC  +  ^HC 


(174) 


The  sign  of  either  or  D^,  will  indicate  tho  amount  of  each  component 
contained  in  the  detected  wave .  The  only,  problem  which  may  arise  will 
be  for  Dl  =  0;  this  will  mean  that  Ir  =  1^  which  may  happen  for  unpolar¬ 
ized  or  for  circularly  polarized  waves.  For  the  type  of  problem  being 
studied  here,  the  U  and  V  parameters  are  always  small  and  circular 
polarization  will  not  result  for  a  transmitted  unpolarized  wave. 

5.  MIE  AND  RAYLEIGH  SCATTERING 


In  general,  the  scattered  components  of  the  electric  field  are 
given  by  [1,4] 


where  the  incident  field  is  a  monochromatic  plane  wave  propagating  in 
the  positive  z  direction  with  components 


e-ikz  +  iwt 


(175) 


and  the  S..  elements  are  scattering  functions  which  depend  on  the 
material  properties  and  orientation  of  the  scatterer. 
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A  transformation  matrix  can  also  be  defined  [9]  for  the  Stokes 
parameters, 


or 


r\ 
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(177) 


(178) 


where  F  and  F  are  4  x  4  matrixes  and  are  given  by  (see  appendix) 


and 


F  = 


\ 


0 

0 


F  = 


/ 


0 

0 


O' , 

°\ 

0  , 

0  } 

21  » 

-D2i  i 

2 1  * 

S2i/ 

M2, 

0  , 

0  , 

0 

0  , 

Mir 

0  , 

0 

0  , 

q  , 

S2i> 

-D  2i 

0  , 

0  , 

D2i » 

S21 

(179) 


(180) 


The  derivation  of  (178)  and  (179)  in  the  appendix  used  the  fact 
that  for  a  spherical  scatterer  A3  =  At*  =  0.  This  is  equivalent  to 
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setting  S3  =  S4  =  0  in  equation  (175).  Therefore,  we  can  write 


l'ES  ) 
l 

1  s2 

0  ^ 

-ikr  +  ikz 
e 

14  1 

ES  i 

0 

S1i 

ikr 

E^ 

r  1 

1 

i  r 

furthermore ,  from  (180);,  we  know  that 

Mi  =  Aj  At  =  |Aj|2  (182) 
M2  =  A2  A$  =  jA2|2  (183) 
S2i  =  l/2(Ai  At  +  At  A2)  (184) 
D2j  =  -  ]■  (Aj  At  -  At  A2)  (185) 


Also,  comparing  (175')  with  (2)  from  the  appendix,  we  obtain 


-ikr  +  ikz 

hi  -  Si — - -  (186) 

ikr 

_ i kr  +  ikz 

A2  =  S2  S - -  (187) 

ikr 

so  that 


|A,I2  |S,|2 

k2  r2 

(188) 

IM2  ■  |sz|2 

k2  r2 

(189) 

s2,  =  -I—  (S,  Sj  +  St  S2) 

2k2  r2 

(190) 

0  = - - —  (Sl  st  -  st  s2) 

2k2  r2 

091) 
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In  general5  Sj  and  S2  will  be  complex: 

Sj  =  sj  ei01  (192) 

S2  =  s2  e1<72  (193) 

\vhere  si  and  s2  are  real  amplitudes  and  cj.  and  o,2  are  real  phases, 
therefore 


i'Aj2  - - — 

(194) 

k2  r2 

S2 

|A2I2  =  — 7 
k2  r2 

S1  S9 

(195) 

S21  =  ‘  „  COS 

k2  r2 

(o 1  -  02) 

(196) 

si  h 

D  =  al — £  sin 
k2  r2 

(a i  -  o2) 

(197) 

Substituting  equations  (188)  through  (191)  into  equation  (180) 
results  in  the  phase  matrix 


1  |s2l2 

0 

0 

0  \ 

j 

1 

0 

IS'll2 

0 

0 

k2  r2 

0 

0 

Re(SjS2) 

-Im(S iS* ) 

\  0 

0 

Im(S1sS) 

Re(S1SJ)  / 

with  S!  =  $!  (e)  and  S2  =  S2  (e)  functions  of  the  polar  angle.  Note 
that,  according  to  (198),  if  the  incident  wave  is  linearly  polarized 
in  the  plane  of  scatter  the  scattered  wave  will  also  be  linearly 
polarized;  this  agrees  with  the  Mie  theory  [7], 
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All  scatt'erers  treated  in  this  study  are  spherical  and  (198) 
will  play  an  important  role  in  the  results.  From  the  theory  of 
Mie,  we  have  that 


and 


Si  (e) 


2n  i  1 

n  Tn  +17 


(cos  e)  +  b  x 
'  '  n  n 


(199) 


JCi 


oo 

c  /.X  _  V  2n  +  1 

^  (0)  "  fa 


*  K 

n  n 


(cos,  e)  +  an  Tn  (cos  e) 


(200) 


The  detailed  definition  of  Tn,  irn,  an,  and  bn  is  postponed  until  the 
next  I’ecti  on . 

The  phase  matrix  for  Rayleigh  scattering  is 


cos*  e 


R=  3/2  ^ 


0 

0 

0 


0  0 

1  0 

0  ’cos  e 

o  0 


(201) 


Equations  (198)  and  (201)  apply  only  when  the  Stokes  parameters 
are  referred  to  the  plane  of  scatter.  In  general,  the  incident 
wave's  Stokes  parameters  will  not  be  referred  to  the  plane  of  scatter, 
therefore  it  wl  1  be  necessary  to  rotate  the  incident  Stokes  parameters 
into  the  plane  of  scatter  [9],  Referring  to  equations  (128)  through 
(131),  it  is  clear  that  the  I  and  V  parameters  are  invarient  to  a 
rotation  of  axes.  If  we  rotate  the  axes  by  an  angle  <{>  in  the 
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clockwise  direction  (figure  10)  then 


Q'  =  a2  cos  2  6  cos  2  (X  -  <f>)  (202) 

and 

U'  =  a2  cos  2  a  sin  2  (x  -  *)  (203) 

Expanding  (202)  and  (203)  We  obtain 

Q'  =  Q  cos  2  <f>  +  U  sin  2  <J>  (204) 

and 

U'  =  -Q  sin  2  <t>  +  U  cos  2  <j>  (205) 

But,  we  also  have  that 

(206) 

V'  ■  V  (207) 

Q'  =  1'  -  I'  =  (I  -  I  )  cos  2  <j»  +  U  sin  2  <j>  (208) 

A  I  Xr  I 

U'  =  -  (I  -  I  )  sin  2  <jj  +  U  cos  2  <f>  ( 209 ) 

x>  r 

and  solving  for  Ir  and  I'  we  have 

x»  r 


If  =  1/2  1^  (1+cos  2  $)  +  1/2  Ir  (1-cos  2  <f>)  +  1/2  U  sin  2  <j» 

=  1^  cos2  <j>  +  I  sin2  <j>  +  1/2  U  sin  2  <j>  (210) 

If  =  1/2  lz  (1-cos  2  *)  +  1/2  Ir  (1+cos  2  <t>)  -  1/2  sin  2  * 

=  I4  S\in2  4  +  If  cos2  $  -  1/2  sin  2  *  (211) 

The  transformation  law  for  a  rotation  of  axes  is  therefore 
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(212) 


/ 


cos2<4 

sin24i 
-sin  2<j> 
\  0 


sin?<{>  1/2  sin,  2<j> 

cos'3$  -1/2  sin  2$ 
sin  2<j>  cos  ?d> 

0  O' 


Prior  to  applying  the  scattering  matrix,  (198)  or  (201 ) ,  the  incident 
Stokes  parameters  must  be  rotated  into  the  plane  of  scatter  according 
to  (212)  and  rotated  back  into  a  reference  plane  if  desired  after 


application  of  (198) 

• 

O 

CO 

o 

Let  us  write 

/  cos2<}> 

sin2<j> 

1/2  sin  2<f> 

°\ 

L  (♦)  = 

|  sin2<j> 

cos2«j> 

-1/2  sin  2$ 

0 

(213) 

1-sin  2ij > 

sin  2$, 

cos  2<}> 

0 

1 

: 

\  O' 

0 , 

0 

1  / 

as  the  transformation  matrix  for  rotation  in  the  clockwise  direction 
then  l  {^4)  is  the  transformation  matrix  for  rotation  in  the  counter¬ 
clockwise  direction. 

The  contribution  to  the  source  function  due  to  scattering  of  a 
pencil  beam  of  radiation  of  solid  angle  dto'  in  the  direction  (e",  <|>') 
is 


for  Rayleigh  scattering  and 


(214) 


(215) 
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for  Mie  scattering  if  I  (e%  $')  is  referred  to  the  directions  parallel 
and  perpendicular  to  the  plane  of  scattering.  In  general,  I  $') 
is  not  referred  to  the  scattering  plane  and  we  must  therefore  rotate 
the  incident  polarization  components  into  the  plane  of  scatter.  Sup¬ 
pose  that  I  ,(e\  $')  is  referred  to  directions  along  the  meridian 
plane  OBA  in  figure  15  and  at  right  angles  to:  it.  We  use  equation  (213) 
to  transform  i>(e'v  $')  into  the  plane  of  scatter,  and  therefore,  the 
contribution  to  the  source  function  is 

R  (cos  9)=W=£y;^i0's  r).  Ifc  (216) 


2 


Pi-.iw.Q  iq  Polarization  Reference  Planes 
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We  can  then  transform  the  Stokes  parameters  back  into  the  reference 
plane  by  once  again  applying  the  transformation  matrix  such  that 

L  (it  -  B2)  R (cos  0)  l(-6i)  I(e',  <T)  37-  (217) 

is  the  scattered  contribution  except  that  now  the  Stcikes  parameters 
are  referred  to  the  plane  OCA. 

If  we  now  assume  the  several  scattering  events  occur  successively, 
that  contribution  due  to/scattering  would  be 

-Ln  -U  ■-  32n*)'!_s  *  -JLr  .(it  1  e2)  R(cos-  0)  L.i.  .(-B.i.).*  •  • 

i  (8-.  r)  (218) 

where  n  such  scattering  events  occur.  If  Mie  scattering  were  taking 
place  rather  than  Rayleigh  scattering  one  would  simply  use- the 
scattering  function  F  rather  than  R  . 

In  this  section  we  have  shown  how  scattering  affects  the  polari¬ 
zation  properties  of  the  incident  intensity  and  have  derived  the 
scattering  matrixes  for  Mie  and  Rayleigh  scattering.  These  ideas 
will  be  used  in  the  modeling  of  scattering  of  electromagnetic 
radiation  in  the  atmosphere. 
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SECTION  IV 
COMPUTER  SIMULATION 


1 .  INTRODUCTION 

A  Monte  Carlo  computer  simulation  of  electromagnetic  scattering  in 
the  atmosphere,  similar  to  those  proposed  by  Kattawar  and  Plass  [16]- ana 
Collins  and  Wells  [2C3,  is  used  to  predict  the  degree  of  depolarization 
of  an  incident  polarized  electromagnetic  wave.  Basically,  the  computer 
program  is  a  simulation  of  the  equation  of  radiative  propagation  in  the 
atmosphere  as  .derived  in  section  II.  The  atmosphere  is  composed  of  a  ran¬ 
dom  distribution  of  particles  of  different  s\izes;  these  which  are  much 
smaller  than  one  wavelength  are  called  Rayleigh  particles  since  they 
follow  the  a"4  scattering  law  originally  derived  by-  Lord  Rayleigh  in  1871; 
while  those  particles  whose  size  is  comparable  to  one  wavelength  or  larger 
are  called  Mie  particles  after  Gustav  Mie  who  first  derived  a  rigorous 
solution  for  the  diffraction  of  a  linearly  polarized  plane  wave  by  a 
sphere.  Since  we  are  principally  concerned  about  polarization  effects, 
absorption  is  neglected. 

The  atmospheric  densities  considered  here  are  selected  such  that  near¬ 
field  effects  may, <be  neglected.  Also,  since  the  particles  which  make  up 
the  atmosphere  are  randomly  distributed,  all  scattering  will  be  assumed 
to  be  incoherent,  and  therefore,  the  Stokes  parameters  can  be  summed  with¬ 
out  regard  to  phase. 

2.  THE  SIMULATION  MODEL 

Consider  the  integral  equation  which  governs  the  radiative  transfer 
of  electromagnetic  energy  in  the  atmosphere 

I  (t  ,m  ,4*)  =  Iq(t,u,<J>)  +  Is(t,u»4>)  (219) 
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where 


Ip(f,p,<i>)  =  1(0, y,*)  e"T/^ 


and 


2tt  +1  T 

I$ (t ,U *4>)  -  ~  J  Jf  F(y,<fr;y',$') 
0  - 1  *0 


I(t,u'  ;<J>' )  e‘(r  “  t)/u  dt  du'  d*! 

The  direct  intensity,  ID(t,y,f),  can  be  calculated  easily.  The  scattered 
intensity  is  much  more  difficult  to  calculate  since  the  unknown  intensity 
also  occurs  under  the  integral.  Equations  of  the  type  (219)  are  known  as 
nonhomogeneous  Fredholm  integral  equations.  In  an  isotropic  and  homoqenous 
medium  the  polarization  parameters  of  the1  direct  intensity  are  only  atten¬ 
uated  and  the  polarization  state  remains  the  same.  The  polarization 
parameters  of  the  scattered  wave  are,  in  general,  modified  during  the 
scattering  process  by  the  scattering  matrix  F(y,<f>;y:  ,<j>' )  (see  equation 
(179)),  where  p  and  <j>  give  the  direction  of  the  incident  wave  while  p' 
and  <j>'  give  the  direction  of  the  scattered  wave, 
a.  Source  Distribution 

Assume  a  source  which  emits  photons  either  isotropically  or  aniso- 
tropically.  The  probability  that  a  photon  is  emitted  in  the  direction 
(e ,4>)  is  proportioral  to  the  magnitude  of  the  source  intensity  pattern  in. 
the  direction  (e,<j>).  Let  the  intensity  pattern  be  denoted  by  the  function 

G(e,<|>)  (220) 

If  we  normalize  the  intensity  function 
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(221) 


2h_7T 

k  f  f  G(e,*) 
0  0 


sin  e  d  e  d <|>  =  1 


where  K  is  the  normalizing  constant,  then  the  probability  that  a  photon  is 
emitted  in  the  interval  (e  +  de,  <|>  +  d<j>)  is 


P(e,<j>)  sin  >6  d  e  d  <j>  =  75 — M  t  -  - 

tTT  7T 

j  j  G(e ,4>)  sin  e  d  e  d  $ 


(222) 


or  upon  integrating  between  fixed  limits 

e2 

j  J  6(e,<j))  sin  0  d  0  d 


P(0l£0<02,(J>i<<J)<_<}i2)  = 


<J>1  01 


2it  tt 

if  G(e,<j>)  sin  0  d  0  d  4. 
0  G 


(223) 


Equation  (223)  gives  the  probability  of  emitting  a  photon  between  the 
angles  (0i,e2)  and  (<j>i,<fca). 

Suppose  the  distribution  is  isotropic,  then 

G(e,$)  =  1  (224) 

and 


($2  -  <h)  (cos  -  cos  02 ) 
P(0  H0102  »^i£i>li>2)  = - j - - 


(225) 


Using  equation  (223)  any  intensity  pattern  may  be  used, 
b.  The  Extinction  Coefficient 

All  scatterers  can  be  characterized  by  their  scattering  function, 
which  describes  the  scattered  intensity  as  a  function  of  aspect  angle.  If 
Iq  is  the  intensity  incident  on  the  scatterer,  then  the  scattered  intensity 
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i 


is  given  by 


I  = 


I0  ^(®»$) 

k2r2 


(226) 


where  F(e,<j>)  is  a  dimensionless  function  which  describes  the  intensity 
distribution  for  all  directions  (e,<f)  and  k  =  2tt/a. 

We  now  define  the  scattering  cross  section,  o$,  by 


=  V  f  He,*) 

k2  J 

ti) 


doj 


(227) 


where  du  =  sin  e  d  e  d  <t>  is  an  element  of  solid  angle,  og  is  an  equivalent 
area  such  that  the  energy  falling  on  og  is  equal  to  the  energy  removed 
from  the  incident  wave.  Notially  there  is  also  an  absorption  cross 
section,  ag,  which  describes  the  amount  of  energy  removed  from  the 
incident  wave  through  absorption.  The  total  extinction  of  the  incident 
wave  is  defined  by  the  sum  of  the  scattering  cross  section  and  absorption 
cross  section. 


o  .  =  a  +  a 
ext  s  a 


(228) 


In  this  paper  we  ignore  oa  and  let 

u 


0  =  °ext  "  °s 


c.  Path  Lengths 

Radiation  intensity  is  attenuated  exponentially  through  a  medium 
according  to 


I(s)  =  I0  e"s/oext(s) 


(229) 


where  Iq  is  the  incident  intensity,  s  is  the  distance  traveled,  and  0  is 
the  extinction  coefficient  of  the  medium.  The  optical  path  length  is 
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defined  as 


(I 

T(h)  =  j  °bxt  (h'}  dh' 


(230) 


where  t(H)  is  also  the  number  of  mean-free-paths  between  0  and  h.  In 
equation  (229)  the  ratio  s/a(s )  is  the  fractional  number  of  meah-free-paths 
traveled  with  s  =  oext(s)  being  equal  to  one  mean-free-path.  If  we  know, 
the  optical  depth  of  the  atmosphere  as  a  function  of  distance  in  mean- 
free-paths,  then  the  number  of  mean-free-paths  traveled  between  collisions 
is:  distributed  exponentially,  i.e.,  the  probability  of  a  collision  between 
p  and  dp,  where  p  =  s/oext(s),  is 


p(p)  dp  =  e"p  dp 


or 


02 

P(p1<p<p2)  =  j  e"p  dp 
P2 


(231) 


(232) 


Equations  (231)  and  (232)  hold  for  0<p<°°  For  0<p<p  we  use 

’  1  liia  a 


D(pi<P£P2) 


P2 

f  e“p  dp 
pi 


(233) 


Equations  (231)  and  (232)  hold  for  an  atmosphere  infinite  in  extent  while 
equation  (233)  holds  for  a  finite  atmosphere 
a.  Rayleigh  Scattering 

The  term  Rayidigh  scattering  is  used  for  scattering  of  electro¬ 
magnetic  waves  from  particles  which  are  small  compared  to  a  wavelength  or 
more  precisely,  particles  which  are  small  compared  to  the  wavelength  in 
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the  particle,  i.e.,  if  Aq  is  the  wavelength  of  the  electromagnetic  wave 
outside  the  particle  then  the  wavelength  inside  the  particle  is 

A, 


m 


(234) 


where  m  is  the  index  <Vf  refraction  (m  =  c/v)  of  the  particle.  If  the 
particle  is  small  compared  to  Aq,  the  wavelength  outside  the  particle 
but  is  comparable  to  the  wavelength,  A,  in  the  .particle -then. we  .will  be 
in  the  resonance  scattering  region.  Resonance  is  associated  with  the 
modes  of  vibration  of  the  particle. 

For  a  particle  fulfilling  the  above  requirements  the  electric 
field  induced  in  the  particle  can  be  assumed  constant  throughout  the 
particle  and  the  induced  dj pole  .moment  is 


P 


a  t 


0 


(235) 


where  a  is  the  polarizability  of  tl.e  particle  and' tQ  is  the  incident 
electric  field.  If 


t  t„  eiu)t 


(236) 


then 


p  =  p0  eiwt  (237) 

The  electric  field  components  of  an  oscillating  dipole  are  given 


Er  =  {{ 


— )  2  cos  8  |p|  ei“(t-R/v) 

R27 

(- —  —  -  — ^  s1„  e  | p I  eia,(t"R/v) 
\R3  R2  R  / 


(238) 


(239) 


;For  distances  large  compared  to  a  wavelength  (R^x)  the  radial  component. 


F  ,  may  be  neglected  and  EA  becomes: 


£  _  k2p  sin  e  e~i kR 


(240) 


:where 

P  -  IPl  e1ut 

Using  (235),  (240)  becomes 


_  k*  «  E„  sin  0 

£ - e 

0  R 


(241) 


We  now  compare  equation  (241)  with  (175)  and  noting  that  for  the  perpen¬ 
dicular  component  o  =  90°  and  for  the  parallel  component  o  =  80c  -  0, 
where  0  is  the  scattering  angle,  we  obtain 


S2  S3^  /cos  0  0 

=  ik3  a 


s.-,  sly 


0  1. 


or 


=  k2  a 


<cos  0  A  -ikr  +  ikz  (  l\ 


0  1/ 


from  which  we  note  that 


(242) 


(243) 


Is  «  |E  r 

i  1  r 


—  cos2  o  jEjl2 
r2  s 


(244) 
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(245) 
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kV 
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2 


2^  cos  <5  cos  0  =  us  cos  0 


i  r 


(246) 


i*  0  =  Js^ai  vs 


Therefore with 
scattering,  becomes 


yipj  sin  6  cos  o  =  r  cos  e  (247) 

1^  «  j^»  U1,  V1'),  the  scattering  matrix,  for  Rayleigh 


R  = 


0  0  0  \ 

1  0  0 

0  cos  0  0 

0  0  cos  0 


(248) 


We  can  now  write  the  matrix  relation  , between  the  incident  and  scattered 
intensity 


Is  =  -^1 R  I1 


(249) 


which  clearly  shows  the  dependence  encountered  with  Rayleigh  scatter¬ 
ing  (k  =  2-n/\). 

e.  Mie  Scattering 

The  exact  solution  to  the  diffraction  of  a  linearly  polarized 
plane  wave  by  a  sphere  was  first  given  by  Gustav  Mie  in  1908  [2].  The 
scattering  of  electromagnetic  waves  by  a  sphere  has  occupied  a  large 
portion  of  the  literature  in  the  past.  Several  forms  of  the  solution 
have  been  given  by  various  authors  [7,  23].  The  term  Mie  scattering  has 
been  used  for  this  problem  in  honor  of  M’e.  in  this  paper  we  will  be 
using  the  solution  given  by  Van  De  Hulst  [14].  The  electric  field 
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components  as  a  result  of  scattering  from  a  sphere  are 

Ee  =  "  ~  S2 ( 9 )  cos  4>  e“lkr+la3t  (250) 

Icr —  ~  - 

and 

%  =  —  S,Ce)  sin  $  e'lkr+ltot  (251) 

y  kr 

where  Si(e)  and  $2(0)  are  the  scattering  functions  defined  by  equation 
(175)  and  are  given  by 


00 


Sl(9)  n(n+1)  j-Vn  <c°'  8>_+  Vn  <cos  e>_ 

(252) 

and 

CO 

S2(o)  '  lj  n(nt])  (Vn  <cos  e>  +  Vn  (cos  6>] 

(253) 

with 

*n  (cos  0)  =  — —  (cos  e) 

sin  e  n 

(254) 

Hn  (cos  0)  =  —  (cos  0) 

"  d0  n 

(255) 

n  ^(y)  ?n(x)  -  m(pn{y)  t^(x) 

(256) 

u  _  "^(y)  *n(x)  -  l«n(y)  ♦'(x) 

"  m^(y)  tn(x)  -  *n(y)  t‘(x) 

(257) 
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In  equations  =(254)  and  (255)  (cos  e),  is  the  associated  Legendre 
polynomial.  In  equations  (256)  and  (257)  *n(x)  and  ?n(x)  are  the  Riccati- 
Bessel  functions  given  by 


tnM  =  x  Vx) 

(258) 

?n(x)  =  x  h^(x) 

(259) 

where  Jn(x)  and  h^(x)  are  spherical  Bessel  functions.  Also  in  equations 
(256)  and  (257)  m  is  the  complex  index  of  refraction  of  the  sphere. 

Equations  (250)  through  (257)  constitute  a  solution  to  the  prob¬ 
lem  of  the  diffraction  of  a  linearly  polarized  plane  wave  by  a  sphere  and 
are  valid  for  an  arbitrary  homogeneous  sphere. 

The  scattering,  matrix  for  Mie  scattering  can  be  obtained  by  sub¬ 
stituting  equations  (252)  and  (253)  into  equation  (198).  Recall  that 
si (e)  is  the  scattering  function  for  perpendicular  incident  polarization 
while  S2(e)  is  the  scattering  function  for  parallel  incident  polarization, 
f.  Description  of  the  Atmosphere 

In  this  study  we  assume  that  the  atmosphere  consists  of  scatter¬ 
ing  particles  oniy.  The  particles  are,  as  is  customary,  assumed  to  be 
spherical.  We  use  the  Mie  scattering  theory  to  determine  the  particle 
scattering  matrix.  Absorption  is  assumed  to  be  negligible  since  this 
phenomenon  does  not  contribute  to  the  depolarization  of  the  incident  wave. 
The  atmosphere  is  made  up  of  many  mie  particles  each  randomly  oriented 
with  respect  to  other  Mie  particles  in  the  atmosphere.  Therefore,  all 
scattering  is  assumed  to  be  incoherent.  The  assumption  of  incoherent 
scattering  allows  us  to  add  the  Stokes  polarization  parameters  resulting 
from  the  various  scattering  events.  Since  the  atmosphere  is  made  up  of 
many  randomly  oriented  Mie  particles,  we  can  describe  the  size  distribution 
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of  the  Mie  particles  and  determine  the  average  scattering  matrix  for  a 
given  size  distribution  and  thus  characterize  the  atmosphere  with  this 
average  value.  Deirmendjian  [24]  has  generalized  the  size  distribution 
function  first  proposed  by  Khrgian  and  Mazin  [25]  in  the  form 

n(r)  =  a  r“  e"brY  (260)' 

where  n(r)  is  the  volume  concentration  at  the  radius  r  and  a,  a,  b,  and  y 
are  positive  constants.  If  N  is  the  total  number  of  particles  per  unit 
volume,  then 

N  =  J  n(r)  dr  =  ^  b-  ~T  (261) 

0 

where  r  is  the  gamma  function. 

Taking  the  derivative .Of. (26QI  with  respect  to  r,  we  obtain 

—  n(r)  =  a  r01”1  e“brY  {«  -  ybrY)  (262) 

dr 


Setting  (262)  equal  to  zero  we  can  solve  for  the  critical  or  mode  radius 
which  is 


Figure  17  is  a  general  plot  of  equation  (260). 


The  generalized  size  distribution  function  given  by  Deirmendjian 
(equation  (260))  reproduces  practical  particle  distributions  as  measured 
by  experimenters.  For  example,  if  a  =  2.373,  a  -  6.,  b  =  1.5,  and  y  =  T 
we  obtain 

n(r)  =  2.373  r6  e"1,5r  cm'3  g'1  (264) 

which  follows  the  size  distribution  measured  by  Durbin  [26]  for  cumulus 
clouds;  for  this  reason  the  distribution  (264)  is  called  "cumulus  type  I" 
by  Singleton  and  Smith  [27],  Another  example  of  equation  (260)  is  obtained 
by  setting  a  =  5.33  x  104,  a  =  1,  b  =  8.944,  and- y  =  1/2  and  we  have 

n(r)  =  5.33  x  104  r  e"8‘944rl/2  (265) 

The  distribution  given  by  equation  (265)  may  be  taken  to  represent 
coastal  conditions  and  is  called  the  "haze  M"  model  by  Deirmendjian. 
Equation  (264)  yields 


N  =  100  particles  cm“3  (266) 

while  equation  (265)  also  gives 

N  =  100  particles  cm"3  (267) 

Equations  (266)  and  (267)  yield,  on  the  average  an  interparticle  distance 
of  approximately 

3 

=  0.215  cm 

or  2150  microns.  Therefore,  for  the  Cumulus  Cloud  and  "Haze  M"  models 
the  far  field  scattering  approximation  is  certainly  valid  for  wavelengths 
up  to,  say,  200g. 
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haze  is 


A  size  distribution  which  may  be  used  to  represent  a  continental 


r  <  r  . 
min 


n(r)  =  0 

n(r)  =  constant  r,,.  <  r  <  r-> 
min  -  -  1 

n(r)  =  Ar"v 


r5  <  r  <  r 
2  -  -  max 


(268) 


where  r  .  =  0.3p,  r2  =  0. 1 y ,  v  =  4,  and  constant  =  103  has  been  shown  by 
Deirmendjian  to  be  representative  of  typical  continental  hazes  and  is 
cabled  the  "Haze  C"  model,  i.e., 

n(r)  =0  r  <  0.03u 

•n(r)  =  103  _0.,.0.3u  £  r  <_  0 . 1  y 

n(r)  =  0.1r“4  O.ly^r^®  (269) 

for  this  distribution  we  also  have 


N  «  100  particles  cm”3 

and  near  field  scattering  will  start  to  become  important  for 

X  >  200y 

Therefore,  it  appears  that  for  wavelengths  info  the  infrared  and  shorter 
near  field  scattering  may  be jgnored  for  the  three  postulated  models. 

Now  that  we  have  obtained  the  three  particle  size  distributions 
we  can  now  integrate  the  Mie  functions  over  the  applicable  size  distribu¬ 
tion  in  order  to  obtain  the  Mie  volume  functions.  Here  we  will  follow 
the  same  technique  used  by  Deirmendjian  [24], 

In  equation  (2?7)  we  defined  the  scattering  cross  section  as 


-  f  F(e,*) 

k2  J 


dw 


(227) 
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where  the  integration  is  carried  out  over  all  solid  angles  w.  If  we 
integrate  the  particle  cross  section  over  ths  size  distribution  we  will 
obtain  the  volume  scattering  function.  Before  carrying  out  the  integra¬ 
tion  let  us  define  the  Mi e  size  parameter 

x  =  ~L  =  kr  (270) 

where  r  is  the  Mie  particle  radius,  k  =  2-rr/x »  and  x  is  the  wavelength. 
Making  this  substitution  into  (260)  we  obtain 

n(x,k)  B  a  j^~J  e" 

The  function  F(e ,4>)  is  also  dependent  on  x  and  we  can  write 

F(e,+)  =  F(e  ,<(>,x )  (272) 

Now  we  can  write,  for  azimuthal! y  symmetric  F(e), 

x2 

2  (e,x)  =  —  f  F(e,x)  n(x,k)  dx  (273) 

s  k2  J 


The  volume  scattering  cross  section  is  obtained  by  integrating  (273)  over 
all  angles,  i.e., 


^  71  ^  ^ 

<^-'vs(m,x)  ~  2vJ  ^^g(0,x)  sin  0 


de 


(274) 


where  m  is  the  index  of  refraction. 

We  can  now  normalize  £s,  i.e., 


Ns(e,x)  =  ^ 


vs 


(275) 
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which  is  proportional  to  the  probability  per -steradian  of  a  photon  of 
wavelength  A  being  scattered  into  an  angle  e. 

For  Mie  scattering  we  have  that 

F(e)  =  Isj2  (276) 

for  incident  perpendicularly  polarized  waves,  and 

F(e)  =  |S2!2  (277) 

for  incident  parallel  polarized  waves  where  the  polarization  is  referred 
to  the  scattering  plane,  i.e.,  either  perpendicular  or  parallel  to  it. 

For  exponential  stratification  of  the  atmosphere  [15]  we  can 
write  the  extinction  coefficient  as 

°ext{h)  =  °0  e“h/H  (278) 

where  oQ  is  the  extinction  coefficient  at  h  =  0  and  H  is  the  scale  alti¬ 
tude,  e.g.,  H  =  0.98  km  for  haze  [15].  Using  equation  (230)  we  obtain 
for  the  optical  thickness 

h 

t(h)  =  J aQ  e'h7H  d  h*  (279) 

0 

Or  upon  solving  equation  (279)  we  have 


which  gives  the  optical  thickness  as  a  function  of  altitude  h. 
g.  A  Monte  Carlo  Simulation 

Monte  Carlo  techniques  may  be  used  in  the  simulation  of  a  random 
process  [12,  13].  Monte  Carlo  techniques  have  been  used  for  several  years 
to  simulate  problems  in  neutron  diffusion.  The  Monte  Carlo  technique 
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serves  in  modeling  a  random  process  when  the  various  probability  distri¬ 
butions  are  known  or  can  at  least  be  approximated  using  physical  reasoning. 
Our  problem  of  determining  the  depolarization  effects  due  to  scattering 
in  the  atmosphere  is  certainly  a  random  process  since  the  distribution 
of  scatterers  in  the  atmosphere  is  random.  Here  we  assume  that  a  plane, 
polarized  wave  propagates  through  the  atmosphere,  the  wave  beinq  made  up 
of  a  large  number  of  photons  whose  average  properties  can  be  characterized 
in  terms  of -wave  theory.  We,  as  originally  stated,  will  also  assume  that 
all  scattering  is  incoherent  so  that  the  effects  of  all  the  scattered 
photons  may  be  added  directly.  The  polarization  properties  of  the  wave 
are  associated^ with  each  photon  history.  The  difficulty  of  assigninq 
the  Stokes  polarization  parameters  to  a  single  photon  can  be  overcome  by 
assuming  that  we  are  dealing  with  groups  of  photons  rather  than  with  a 
single  photon. 

We  now  know  that  the  probability  of  a  photon  being  emitted  into  an 
angle  (e,<}>)  is  given  by  equation  (223),  that  the  number  of  mean-free-paths 
traveled  between -collisions  is  given  by  the  distribution  (233).  Further¬ 
more,  if  Rayleigh  scattering  occurs  the  Rayleigh  scattering  matrix  given 
by  (250)  can  be  used  to  operate  on  the  Stokes  vector  as  given  by  equation 
(249).  Similarly,  if  Mie  scattering  occurs  the  Mi e  scattering  matrix 
given  by  equations  (198),  (252),  and  (253)  can  be  used  to  operate  on  the 
Stokes  vector  as  in  the  case  for  Rayleigi.  scattering.  Since  the  Rayleigh 
and  Mie  scattering  matrixes  are  given  in  terms  of  the  scattering  angle, 
the  polarization  components  of  the  incident  photon  must  be  rotated  as 
dictated  by  equation  (250)  parallel  and  perpendicular  to  the  scattering 
plane.  Similarly,  after  the  scattering  matrix  is  applied  to  the  Stokes 
vector,  the  Stokes  components  will  be  rotated  back  to  a  reference  plane. 

The  reference  plane  to  be  used  here  is  the  same  as  that  used  by  Collins 


and  Kells  [17]  and  is  the  vertical  plane  which  also  contains  the  propaga¬ 
tion  vector.  Therefore,  the  parallel  component  at  a  detector  is  that 
component  in  the  reference  plane  and  is  orthogonal  to  the  propagation 
vector,  while  the  perpendicular  component  will  be  perpendicular  to  both 
the  reference  plane  and  the  propagation  vector. 

Again  following  the  same  method  as  that  used  by  Collins  and  Wells, 
we  determine  the  effect  of  multiple  collisions  by  estimating  the  intensity 
at  the  detector  after  each  collision  and  comparing  the  intensity  at  the 
detector  due  to  single  or  multiple  collisions.  The  estimate  of  the 
intensity  at  the  detector  as  a  function  of  collision  number  will  be 
obtained  by  calculating  the  scattering  angle,  toward  the  detector  and 
weighting  the  Stokes  vector  with  the  applicable  scattering  matrix  for 
that  particular  scattering  angle.  The  particular  photon  history  is  con¬ 
tinued  by  selecting  a  scattering!  angle  from  a  probability  density  given 
by  the  normalized  Rayleigh  or  Mie  scattering  functions.  After  selection 
of  the  scattering  angle  the  Stokes  vector  is  then  weighted  by  the  applic¬ 
able  scattering  matrix  for  that  scattering  angle. 

Random  numbers  are  obtained  by  generating  a  pseudorandom  sequence. 
Several  methods  for  generating  pseudorandom  numbers  are  available  [12]. 

In  tfrs  simulation  we  will  use  a  pseudorandom  sequence  which  is  generated 
by  the  recurrence  relation 

yi  =  ayi  -  1  (modulo  n)  (281) 

The  notation  used  in  (281)  means  that  y.  is  the  reminder  when  ay^  -  1  is 
divided  by  n.  In  our  case  we  use  n  =  248  and  the  period  is  n/4,  hence  we  ’ 
are  able  to  generate  in  the  order  of  1013  numbers  before  repeat-. ng  the 
sequence.  We  will  never  require  more  than,  say,  104  photon  histories 
each  requiring  approximately  102  random  numbers,  thereto 'e  we  will  never 
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need  more  than  10b  random  numbers  per  problem.  The  pseudorandom  sequence 
given  by  (281)  has  been  found  to  be  uniformly  distributed  and  therefore 
a  transformation  from  the  uniform  distribution  to  the  desired  distribution 
is  needed. 

A  theorem  by  Barkovskii  and  Smirnov  [13]  states  that  if  the 
random  variable  z  has  a  probability  density  function  f(x),  then  the 
distribution  of  the  random  triable 

n  =  f  f (x)  dx  (282) 

0 


is  uniform  in  the  interval  (0,1).  Consider  the  arbitrary  cumulative 
distribution  function  F(0  as  given  in  figure  18. 


Let  n  =  FU);  we  now  want  to  find  out  how  n  is  distributed.  Since  we 
know  that  the  aridom  variables  z  and  n  are  functionally  related  through 
the  function  F  we  can  write 


But  we  know  that 


p(n) 


P(6) 

f'u; 


(283) 


pU)  =  F'U) 


(284) 


therefore 


p(o)  =1  0  <  n  £  1 

p(n)  =  0  otherwise,  which  Droves  that  n  is  uniformly  distributed.  On  the 
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basis  of  equation  (282)  we  may  therefore  use  the  following  rule.  In  order 
to  obtain  a  random  number  belonging  to  the  set  of  numbers  |s.j|  withthe 
density  function  f(x)  solve  the  following  equation  for 


/ fW 


dx  =  R. 


(285) 


For  example  consider  the  distribution  function  given  by  equation  (232) 


p(p)  dp  =  e"pdp,  0  <_  P  <_ 


(232) 


which  gives  the  probability  of  a  photon  traversing  dp  mean-free-paths 
between  collisions.  Assume  we  generate  a  random  number  R.  from  a  rectan¬ 
gular  distribution  in  the  interval  (0,1),  then 


Ri  ■  J  e‘P  d> 


=  1  -  e’pi 


Solving  for  p.. ,  we  obtain 


p«  =  -an 


(,  -  Rt) 


(286) 


(287) 


In  this  manner  we  can  very  easily  generate  the  random  numbers  p.  for  the 
desired  distribution  from  the  random  numbers  R.  which  are  uniformly 


distributed. 


Let  us  now  summarize  the  methods  being  used  to  model  our  problem. 
Consider  a  plane-parallel  atmosphere  with  exponential  particle  stratifica¬ 
tion  A  source  of  electromagnetic  energy  is  located  at  an  arbitrary  alti¬ 
tude  above  the  ground.  The  electromagnetic  energy  being  radiated  from 
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need  more  than  106  random  numbers  per  problem.  The  pseudorandom  sequence 
given  by  (281)  has  been  found  to  be  uniformly  distributed  and  therefore 
a  transformation  from  the  uniform  distribution  to  the  desired  distribution 
is  needed. 

A  theorem  by  Barkovskii  and  Smirnov  [13]  states  that  if  the 
random  variable  5  has  a  probability  density  function  f(x),  then  the 
distribution  of  the  random  variable 

n  =  f  f (x)  dx  (282) 

0 

is  uniform  in  the  interval  (0,1).  Consider  the  arbitrary  cumulative 
distribution  function  F(c)  as  given  in  figure  18. 


Figure  1,8.  Arbitrary  Cumulative  Distribution  Function 
Let  n  =  F'U);  we  now  want  to  find  out  how  n  is  distributed.  Since  we 
know  that  the  random  variables  5  and  n  are  functionally  related  through 
the  function  F  we  can  write 

P«n)  -  (283) 

F’(5> 

But  we  know  that 

p(€)-r*(0  •  /  (284) 

therefore 

p(n)  =  1  0<n<l 

p(n)  a  0  otherwise,  which  Droves  that  n  is  uniformly  distributed.  On  the 
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the  source  will  be  modeled  by  conside:  ing  that  the  intensity  is  due  to 
the  conglomerate  effects  of  wav^  packets  made  up  of  photons.  Each  wave 
packet  is  assumed  to  be  characterized  by  the  Stokes  polarization  vector. 
The  wave  packets  can  be  thought  of  as  bullets  which  when  she'  into  the 
atmosphere  collide  and  are  reflected  by  the  scattering  particles  according 
to  prescribed  probabilities. 

The  probability  th  t  a  photon,  (wave  packet)  is  emitted  by  the 
source  in  the  interval  of  solid  angle  (di<e<62,<P}l<Pl<P2)  is 


4>2  02 


//G(e- 


<f>)  sin  e  d  e  d  <f» 


f>(0ll0l02»  <J>  ll'f’lfo )  = 


<h  0i 


2  7T  7T 

j  j G(e»<{>) 
0  0 


(288) 


sin  0  d  e  d  .j) 


where  G(e^)  is  the  source  intensity  distribution.  Therefore  once  a 
specific  function  G ( e )  is  selected  we  can  generate  random  numbers 
representing  a  random  angle  over  the  prescribed  distribution  through  the 
use  of  equation  (285). 

The  probability  that  the  photon  traverses  the  distance  between 
pj  and  p2  between  collisions  or  from  the  source  to  the  first  collision  is 
given  by 


P(pi<p<p2) 


(289) 


where  p  is  the  number  of  mean-free-paths  traveled.  We  will  use  (287)  to 
generate  the  random  numbers. 

We  can  see  from  figure  19  that  the  exponential  distribution  for 
optical  pathslengths  holds  for  the  interval  0<p£». 
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Figure  19.  Exponent*  ■  *  l.>  Distribution 

If  we  desire  to  limit  the  optical  to  the  interval  0<p<pm3U  as  in  the 
case  of  a  very  dense  c'/oud  iwith^  eligible  scattering  outside  the  cloud 
or  for  a  finite  atmosphere,  a  truncated  exponential  distribution,  as 
depicted  in  figure  19  would  be  sed.  The  truncated  exponential  random 
numbers  may  be  found  from 

i 


dp 


or,  upon  solving  for  p. 


e"pmax 


(290) 


(29T) 


As  we  have  mentioned  previously,  we  will  estimate  the  intensity 
at  a  detector  after  each  collision.  We  do  this  by  calculating  the 
scattering  angle  toward  the  detector v  rotating  the  Stokes  vector  into  the 
scattering  plane,  multiplying  the  scattering  matrix  with  the  rotated 
stokes  vector,  rotating  the  Stokes  vector  back  to  the  reference  plane, 

*nd  attenuating  the  resultant  Stokes  vector  by  the  attenuation  function 


e"pD 

D2 


where  p^  is  the  optical  distance  from  the  collision  point  tu  the  detector 
and  D  is  the  gsomtric  distance  from  the  collision  point  to  the  detector. 
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Figure  20  shows  the  geometry  used  for  calculating  the  intensity 
at  a  detector  resulting  from  a  collision.  In  figure  20*  r  is  the  hori¬ 
zontal  radial  distance  from  the  origin  to  the  collision  point,  rD  is  the 
horizontal  radial  distance  from  the  origin  to  the  detector,  o  is  the 
scattering  angle  toward  the  detector,  <j>r  is  the  Stokes  rotation  angle, 
hD  is  the  height  of  the  detect.dr,  and  hc  is  the  height  of  the  collision 
point.  Using  the  spherical  triangle  age  we  can  compute  the  rotation 
angle  <£r  from 


cos  <f>r 


COS  a  -  COS  0  COS  g 

sin  e  sin  g 


(292) 


sin  $r 


_  sin  a  sin  y 
sin  o 


(293) 


where  y  is  the  angle  between  the  radial  lines  to  the  collision  point  and 
to  the  detector  as  shown  in  figure  20.  The  cosine  of  the  scattering  angle 
may  be  computed  by  taking  the  inner  product  between  the  unit  vector  in  the 
direction  of  the  photon  prior  to  the  collision  with  the  unit  vector  in  the 
direction  from  the  collision  point  to  the  detector.  The  unit  vector  in 
the  direction  of  the  photon  prior  to  the  scattering  event  is 

£p  =  sin  g  cos  6  ax  +  sin  g  sin  6  ay+  cos  g  az  (294) 

where  the  angle  5  is  shown  in  figure  20.  The  unit  vector  from  the  colli¬ 
sion  point  to  the  detector  is 

aD  =  -  (rD  cos  y  -  r)  ax  +  rp  sin  y  ay  +  (hp  -  h)  az  (295). 


Taking  the  inrsr  product  of  equations  (294)  and  (295)  we  obtain 
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Figure  20,  Geometry  for  Calculation  of  Intensity  at  a  Detector 


(296) 


D 


After  having  made  an  estimate  of  the  intensity  at  the  detector 
after  a  collision,  a  random  scattering  angle  is  selected.  The  scattering 
event  may  be  either  Rayleigh  scattering  or  Mie  scattering.  The  specific 
type  of  scattering  which  occurs  for  a  particular  collision  is  determined 
by  comparing  a  random  number  Selected  from  a  uniform  distribution  with 
the  ratio  of  Rayleigh  to  Rayleigh  plus  Mie  scattering  coefficients  for 
that  altitude  [28].  Tf  the  random  number  is  less  than  or  equal  to  the 
ratio  Rayleigh  scattering  is  selected.  Otherwise  we  select  Mie  scattering. 

For  Rayleigh  scattering  the  scattering  angle  is  selected  from  the 
probability  density  function  [9,  20] 


P(cos  0)  s  —  ('1  +  cos2  0) 

<(297) 

8 

Collins  and  Wells  have  shown  that  (297)  can  be  written  as 

the  sum  of  two 

probability  density  functions 

P(cos  0)  =  Pi(co$  0)  +  P2 (cos  0) 

(298) 

with 

Pi (cos  0)  -  2.  s1"  6  d  9 

(299) 

4  2 

and 

P2{C0S  oMcOs2  0Sln9d8 

(300) 

4  2 

The  density  function  Pi (cos  0)  is  chosen  to  represent  P(cos  0)  half  of  the 

time  and  P2(cos  o)  the  remainder  of  the  time.  If  Pi(cos  0)  is  chosen, 
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then  cos  0  may  be  determined  by  solving 


..  I  sin  e  d  e 

1  J.  ? 


(301 ) 


for  cos  Q- ,  or 


cos  0.  =  1  -  2  R.j 


(302) 


If  P 2 (cos  0)  is  chosen,  cos  is  selected  in  the  same  manner  as  (302) 
but  cos2  0^  is  compared  to  a  second  random  number  and  is  accepted  if  it 
is  greater  than  the  second  random  number,  Collins  and  Wells  have  shown 
that  the  determination  of  cos  0  compares  very  favorably  with  the  actual 
distribution  for  Rayleigh  scattering. 

If  Mie  scattering  is  selected,  the  cumulative  probability  of 
scattering  into  ah  angle  less  than  or  equal  to 0^  is  given  by  the  normal 
ized  differential  scattering  cross  section  of  equation  (273): 


>(0<O<Oi)  = 


(0, A )  sin  0  d  0 


(303) 


(o,a)  sin  o  d  0 

0 


The  azimuthal  scattering  angle  which  determines  the  orientation 
of  the  scattering  plane  is  selected  from  the  uniform  distribution 


p(<f>)  =  --  (0<^<2ir) 
2tt 


(304) 


for  both  Rayleigh  and  Mie  scattering. 

Now  that  the  scattering  polar  and  azimuthal  anqles  are  known  we 
my  compute  the  direction  of  the  photon  after  being  scattired.  The  cosine 
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of  the  angle  between  the  photon's  direction  after  scatter  and  the  vertical 
axis  is 


cos  a'  =  cos  81  cos  0  +  sin  8'  sin  0  cos  <jy  (305) 


where  the  angles  are  defined  as  shown  in  figure  21. 

The  intensity  received  at  each  detector  simulated  will  be  recorded 
as  a  function  of  polar  and  azimuthal  angles.  The  polar  axis  for  each 
detector  is  defined  to  be  the  axis  joining  the  source  and  detector.  The 
x  and  y  axes  lie  in  a  plane  normal  to  the  polar  axis.  The  x  axis  is 
contained  in  the  vertical  plane  containing  the  source  and  receiver  points. 
The  positive  x  axis  points  toward  the  positive  vertical  direction  but  is 
in  general,  inclined  to  the  true  vertical. 

If  reflections  from  the  ground  occur,  we  assume  that  the  reflected 
light  is  diffuse  and  therefore  completely  depolarized.  The  reflected 
polarization  parameters  are  given  by 


I 


I 

r 


(307) 


u*  =  0 


(308) 


V  =  0  .  (309) 

where  the  prime  denotes  the  reflected  quantities  and  a  is  the  ground 
alb'edo. 
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SECTION  V 


DEPOLARIZATION  BY  CLOUDS:  SOME  EXAMPLES 
1-../  INTRODUCTION 

Ip  the  previous  sections  we  have  discussed  the  various  facets  of 
the -problem  of  determining  the  effect  of  atmospheric  scattering  on  the 
polarization  properties  of  an  electromagnetic  wave.  In  section  II  we 
derived  the  equation  of  radiative  transfer  for  a  plane-parallel  atmos¬ 
phere.  There  it  was  found  that  the  intensity  at  any  po/nt  may  be 
'found  by  solving  a  complex  set  of  integral  equations.  The  equation  of 
radiative  transfer  has  been  solved  only  for  the  cases  of  isotropic 
scattering  and  Rayleigh  scattering  [9] .  The  more  general  problem  of 
anisotropic  scattering,  has  also  been  discussed  by  Evans,  Chu,  and 
Churchill  [11].  The  scattering  function  for  single  scattering  may  be 
expanded  in  a  series  of  Legendre  polynomials  of  the  form 

N 

f ( cos  0)  =  pi  (cos  0)  (310) 

i=0 

where  the  a.  are  the  series  coefficients.  This  expansion  is  of  the 
same  form  as  the  exact  solution  for  a  sphere  (see  equations  (252)  and 
(253)).  The  scattering  functions  selected  by  Evans  et  al.  are  highly 
peaked  in  the  forward  direction  which  indicates  that  for  that  particu¬ 
lar  case  near  field  effects  would  be  more  pronounced  in  the  forward 
scattering  direction.  The  shape  of  the  scattering  function  will 
depend  on  the  size  parameter  x  =  where  r  is  the  radius  of  the 
spherical  particle  and  X  is  the  wavelength;  it  will  also  depend  on  the 
sphere's  complex  index  of  refraction. 
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2.  PARTICLE  SCATTERING  FUNCTIONS  FOR  THE  CUMULUS  CLOUD  DISTRIBUTION 
In  section  V  we  derived  the  normalized  scattering  function  for  a 
given  particle  size  distribution;  in  this  subsection  we  will  show  some 
results  for  the  "Cumulus-  Cloud"  distribution.  This  distribution  is 

n(r)  =  2.373  r';  e  1>5r  cm"3  u"1  (  11 ) 

as  given  by  Deirmendjian,  In  equation  (311)  the  particle  density  is 

N  =  ,oo  PArMcJ.es 

cm 3 

and  the  mode  radius  is  rm  =  4j.  ,  Equation  (311)  is  plotted  in  figure 
22.  In  figures  23  through  27  are  shown  the  scattering  function  for 
individual  particles  taken  from  the  Cumulus  Cloud  distribution  and 
perpendicular  incident  polarization,  with  X  =  b.5u  and  i^e)  =  ( Sj  (e)  | 2 . 
The  size  parameters  13.86,  35.87,  72.31',  123.19,  and  188.50  correspond 
to  spheres  of  radius  1,1,  2.9,  5.8,  9.8,  and  15. On,  respectively. 

Notice  the  increase  in  the  forward  peak  with  size  parameters. 

We  can  normalize  equation  (311)  as  follows: 

p(r)  -  -."td--  -  ^  (312) 

j  n(r)  dr 
0 

or  upon  substituting  for  n(r)  from  equation  (311)  and  using  the  value 
of  N  we  have 


p(r)  =  0.02373  r6  e~1*5r  (313) 

which  is  the  probability  of  finding  a  particle  in  the  radius  interval 
(r,  r  +  dr).  The  expected  value  of  r  can  be  found  from 
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E[r]  =  J  r  p(r)  dr  (314) 

0 

or 

oo 

E[r]  =  f  0.02373  r  e~1,5r  dr 
J: 

o 

=  0.02373  (l.i>)“8  r(8) 

=  4.66m  (315) 

The  variance  is 

VarO]  =[e  (r  -  E[r])2J 
or 

Var[r]  .«  E[r2]  -  (E[r])2  (316) 

,le  can  obtain  the  second  term  on  the  RKS  from  (315)  and  the  first  term 
from 

co 

E[r2]  =  J  r2  p(r)  dr  (317) 

0 

or 

CO 

E[r2]  =  J  0.02373  r8  e"1*5r  dr 
0 

=  0.02373  (1.5)“9  r  (9) 

=  24.9m2  (318) 

so  that  upon  using  the  result  (315),  we  obtain 
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Var[r]  =  24.9  -  (4.66)2 
=  3.2p2 


(319) 


We  can  now  obtain  the  standard  deviation 

a  =  x/VarCr]  =  1.78m  (320) 

Comparing  the  expected  value  of  r,  E[r],  with  the  mode  radius,  r  ,  we 
find  that  since  these  values  do  not  coincide  and  the  distribution  is 
not  symmetrical  the  particles  of  radius  greater  than  E[r]  will  affect 
the  averaging  of  the  scattering  function  over  the  size  distribution  to 
a  greater  extent  than  those  below  E[r].  Comparing  figures  23  through 
27  we  can  see  that  the  scattering  functions  given  by  figures  24  through 
26  will  affect  the  averaging  of  the  scattering  function  over  the  size 
distribution  to  a  greater  extent  than  the  scattering  functions  given 
by  figure  23  or  figure  27.  The  scattering  function  given  by  figure  25 
will  enter  into  the  averaging  quite  heavily  since  it  is  within  la  of 

V 

E[r].  The  fact  that  particles  of  radius  greater  than  rm  will  greatly 
affect  the  normalized  scattering  function  can  be  seen  by  comparing 
figures  25  and  26  with  figure  28.  Figure  28  is  the  normalized  scat¬ 
tering  function  for  perpendicular  incident  polarization  obtained  from 
equation  (275).  The  amplitudes  of  the  normalized  function  and  indi¬ 
vidual  scattering  function  cannot  be  compared  since  the  individual 
particle  scattering  functions  have  not  been  normalized.  The  normalized 
scattering  function  for  parallel  incident  polarization  is  shown  in 
figure  29.  There  is  not  much  difference  between  the  normalized  scat¬ 
tering  functions  for  perpendicular  and  parallel  polarization  except  in 
the  scattering  interval  from  90°  to  170°.  The  end  points  at  scattering 
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angles  of  0°  and  180°  of  the  normalized  scattering  functions  for 
parallel  and  perpendicular  incident  polarization  are  equal  as  required, 
i.e.,  the  energy  scattered  in  the  forward  or  back  directions  are  not 
sensitive  to  polarization. 

If  we  average  the  normalized  scattering  functions  for  perpendicular 
and  parallel  incident  polarization  as  follows: 

I ( e) | 2  +  | S  0 ( e ) | 2  ij(0)  +  l?(e) 

— - - - -  - - 2“^”  (S&l) 

we  will  have  ,a  relation  which  is  proportional  to  the  probability  of 
scattering  in  any  direction  and- we  can  derive  the  cumulative  scatter¬ 
ing  distribution.  This  vividly  shows  (figure  30)  the  preference  to 
forward  scattering  for  the  cumulus  cloud  distribution  at  a  wavelength 
of  0.5p. 

If  we  now  reduce  the  frequency  (increase  the  wavelength)  and  con¬ 
tinue  to  use  the  same  size  distribution  given  for  cumulus  clouds  and 
% 

consider  particles  of  the  same  radius  as  before  we  can  get  a  "feel" 
for  the  effect  of  wavelength  on  the  scattering  function.  Figures  31 
through  33  give  the  scattering  functions  for  particle  radii  of  1 J, 

2.9,  5.8,  9.8,  and  15. 0y  as  before.  These  are  given  for  perpendicular 
incident  polarization  and  a  wavelength  of  4.0 u.  Figures  36  and  37 
give  the  normalized  scattering  functions  for  perpendicular  and  parallel 
incident  polarizations,  respectively.  Figure  38  is  the  cumulative 
scattering  function.  Again  we  see,  upon  comparing  figures  33  and  34 
with  figure  36,  that  radii  greater  than  rm  affect  the  normalized  scat¬ 
tering  function  to  a  great  extent.  We  can  also  compare  figures  30  and 
38  and  see  hhat  the  probability  of  scattering  in  the  forward  direction 
when  compared  to  side  and  back  scattering  has  been  reduced  somewhat. 
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We  continue  to  increase  the  wavelength  to  6. Op  and  then  to  10. On. 
Figures  39  through  43  and  figures  47  through  51  are  the  scattering 
functions  for  perpendicular  incident  polarization  for  wave  lengths  of 
6.0  and  lO.Oy,  respectively  and  given  at  the  same  radii  as  in  figures 
23  through  27  and  figures  31  through  35.  The  widening  of  the  scat¬ 
tering  functions  with  increasing  wavelength  can  be  seen  by  comparing 
figure  30,  the  cumulative  probability  function  for  X  =  0.5m  with 
figure  46,  the  cumulative  probability  function  for  X  =  6,0m  and  fig¬ 
ure  54,  the  cumulative  probability  function  for  X  =  10.0m.  For 
example,  in  figure  30,  the  probability  of  scattering- through  an  angle 
between  0°  and  21°  is  0.5,  the  same  probability  as  obtained  for  the 
scattering  interval  from  4°  to  21°  in  figure  54.  This  dependence  of 
scattering  angle  on  wavelength  implies  that  for  a  given  particle  size 
distribution  that  multiple  scattering  will  increase  with -wavelength, 
i.e.,  the  probability  of  scattering  over  a  greater  range  of  angles  is 
greater  for  large  wavelength.  Figures  44  and  45  are  the  normalized 
scattering  functions  at  X  =  6.0m  while  figures  52  and  53  are  the 
normalized  scattering  functions  for  X  =  10.0m. 

The  extinction  coefficient  for  Mie  scattering  for  the  cumuius  cloud 
distribution  is  given  in  table  1  for  each  wavelength  under  considera¬ 
tion.  The  extinction  coefficient  for  Mie  scattering  is  given  by  [14] 


°ext  <x-r'n'>=lvE<2"  +  1>Re(an  +  bn) 
n=l 


(322) 


where  the  coefficients  an  and  bn  are  given  by  equations  (256)  and  (257). 


We  can  define  the  extinction  efficiency  factor  by 


Wx’m> = 


°ext 


nr" 


(323) 
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Figure  23.  Particle  Scattering  Function  for  X  =  0.5u  and  X  =  13. 86, 

Cumulus  Cloud  Distribution 
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Figure  24.  Particle  Scattering  Function  for  X  =  0.5y  and  X  =  35.07 
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Figure  25.  Particle  Scattering  Function  for  X  =  0.5u  and  X  =  72.31 
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Figure  27.  Particle  Scattering  Function  for  X  =  0.5p  and  X 
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Figure  28.  Volume  Scattering  Function  I x (e )  for 
Cumulus  Cloud  Distribution  at  X  =  0.5y 
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f  ’  Figure  30.  Cumulative  Distribution  for  Cumulus  Cloud 

i  Distribution  at  X  =  0.5;t 
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Figure  32.  Particle  Scattering  Function  for  X  =  4, On  and  X  =  4.48 
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Figure  36.  Volume  Scattering  Function  Ij(e)  for 
Cumulus  Cloud  Distribution  at  \  =  4.0m 
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Figure  37.  Volume  Scattering  Function  I2(0)  for 
"Cumulus  Cloud  Distribution  at  X  =  4. Op 
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CUMULATIVE  DISTRIBUTION 


Figure  41.  Particle  Scattering  Function  for  \  3  6.0u 
and  x  =  6.03,  Cumulus  Cloud  Distribution 
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Figure  44.  Volume  Scattering  Function  1 1 (e )  for  Cumulus  Cloud 
Distribution  at  x  =  6. Op 


SCATTERING  FUNCTION  I2(e) 


Figure  45.  Volume  Scattering  Function  I2(e)  for  Cumulus  Cloud 

Distribution  at  \  =  6.0y 
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CUMULUS  CLOUD  DISTRIBUTION- 
WAVELENGTH  =  6.00  MICRONS 


SIZE  PARAMETER  =  .69 
WAVELENGTH  =  10.00  MICRONS 
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Figure  47.  Particle  Scattering  Function  for  \  -  10. Oy 
and  x  =  0.69,  Cumulus  Cloud  Distribution 


Figure  49.  Particle  Scattering  Function  for  X  =  10. Op  and  X  =  3.6 
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Figure  50.  Particle  Scattering  Function  for  X  =  10. Op  and  X  =  6.16, 

Cumulus  Cloud  Distribution 
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Figure  52.  Volume  Scattering  Function  I x (e)  for  Cumulus  Cloud 

Distribution  at  \  =  10. Ou 


121 


SCATTERING  FUNCTION  I~(0) 


Figure  53.  Volume  Scattering  Function  12(0)  for  Cumulus  Cloud 
Distribution  at  x  s  10, Oy 
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Figure  54.  Cumulative  Distribution  for  Cumulus  Cloud 
Discribution  at  \  =  10. On 


Table  1 


EXTINCTION  COEFFICIENTS  FOR  THE  CUMULUS  CLOUD  DISTRIBUTION 


X(n) 

2ext 

0.5 

1.657  x 

10'4 

1.0 

1.712  x 

10-4 

2.0 

1.812  x 

10‘4 

4.0 

2.179  x 

10"4 

6.0 

2.537  x 

10“4 

8.0 

2.303  x 

10“4 

10.0 

1.865  x 

10“4 

where  r  is  the  particle  radius.  Then  by  integrating  over  the  size 
distribution  [24]  we  obtain  the  volume  extinction  cross  section  as 

r 


^ext  (A 


•">*  / 


7T r2  Q 


ext 


(H n(r) 


dr 


(324) 


Since  we  are  not  considering  complex  indexes  of  refraction  in  this 
study  the  extinction  cross  section  is  equal  to  the  scattering  cross 
section.  We  also  notice  on  comparing  the  extinction  coefficients  in 
table  1  for  each  wavelength  that  they  are  almost  constant  as  a  function 
of  wavelength. 

3.  PARTICLE  SCATTERING  FUNCTIONS  FOR  HAZE  M  DISTRIBUTION 
We  now  consider  the  particle  size  distribution  (figure  55) 

-  0  ,  qtiL 

n(r)  =  5.33  x  104  r  e  (325) 
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which  is  typical  of  a  Maritime  Haze  [24].  As  before  we 
distribution 


P(r) 


Mr) 

N 


where  again  N  =  TOO 
and 


particles }  and 
cm3 


p(r)  =  533  r  e"8’944^ 


The  expected  value  is 


oo 

E[r]=f  533  r2  e”8,944t/^  dr 
0 

=  533  (8.944)"6  r(6) 

=  0 . 1 25y 


The  mode  radius  is 


m 


0.05y 


The  variance  is 

Var[r]  =  E[r2]  -  (E  [r])2 


wi  th 


and 


oo 

E[r2]  =  j  533  r3  e“8‘944^'dr 
0 

=  533  (8.944)” 8  r(8) 

=  0.07y2 


Var[r]  =  0.07  -  (0.125)2 
=  0.05y2 


normalize  the 

(326) 

(327) 

(328) 

(329) 

(330) 

(331) 
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The  standard  deviation  is  then 

aH  =  ^  Var[r]  =  0.22u  (332) 

Comparing  the  mode  radius  with  the  expected  value  we  can  see  that 
particles  with  radii  greater  than  the  mode  radius  will  again  contribute 
heavily  to  the  normalized  scattering  function. 

We  can  compare  the  ratios  of  the  standard  deviation  to  the  expected 
value  for  the  Cumulus  Cloud  distribution  and  for  the  Haze  M  distribu¬ 
tion.  We  have  for  the  cumulus  distribution 


°c  _  1.78 
E[rJ  "  4.66  “ 


0.38 


while  for  the  Haze  M  distribution 


(333) 


qH  _  0.22  _  ,  7, 

UrT  '  0.125  “  1,76  (334) 

from  which  we  can  see  that  the  relative  spread  is  larger  for  the  Haze 
M  distribution  and  we  conclude  that  a  wider  range  of  particle  sizes 
will  contribute  to  a  greater  extent  for  the  Haze  M  model  than  for  the 
Cumulus  model. 

Figures  56  through  60  show  the  scattering  functions  for  perpendicu¬ 
lar  polarization  and  X  =  0.5m  for  individual  particles  of  radius  0.09, 
0.33,  0.73,  1.29,  and  2.00m  with  size  parameters  of  1.06,  4.09,  9.12, 
16.13,  and  25.13,  respectively.  Here  again  we  see  an  increase  in  for¬ 
ward  scattering  with  increasing  size  parameter.  Figures  61  and  62  are 
the  normalized  scattering  functions  for  perpendicular  and  parallel 
incident  polarizations,  respectively.  In  comparing  figures  55  through 
60  with  figure  61  we  can  see  that  the  larger  particles  contribute  to 
the  normalized  scattering  function  more  than  do  the  smaller  particles. 
The  cumulative  probability  distribution  is  shown  in  figure  63, 
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Comparing  figures  30  and  63  we  see  that  the  scattering  angles  are 
distributed  over  a  slightly  wider  range  for  the  Haze  M  model  than  for 
the  Cumulus  model,  ->.6.,  the  probability  of  scattering  between  the 
angles  0°  to  17°  is  0.5  for  Haze  M  while  for  Cumulus  for  the  same 
probability  the  interval  is  between  0°  and  4°.  Therefore,  for  a  given 
wavelength  and  density  we  would  expect  more  multiple  scattering  for 
Haze  M  than- for  Cumulus. 

As  we  increase  the  wavelength  to  2. Op  the  particle  scattering 
functions  increase  significantly  as  shdwn  in  figures  64  through  68 
which  results  in  a  fairly  wide  normalized  scattering  function  as  shown 
in  figures  69  and  70.  The  cumulative  probability  function  is  shown  in 
figure  71 . 

The  particle  scattering  function  continues  to  widen  (figures  72 
through  76)  for  \  =  6. Op,  with  figure  72  showing  an  almost  isotropic 
scattering  function  for  x  =  0.09  and  a  radius  of  0.09p.  The  normalized 
scattering  functions  in  figures  77  and  78  show  this  increase  in  width 
as  does  the  cumulative  probability  function  in  figure  79. 

At  a  wavelength  of  10. Op  the  normalized  scattering  function  for 
perpendicular  incident  polarization  (figure  85)  is  almost  isotropic 
and  the  normalized  scattering  function  for  parallel  incident  polariza¬ 
tion  is  almost  equal  in  the  forward  and  back  directions  with  a  null 
at  90°.  This  type  of  scattering  function  is  typical  of  Rayleigh  scat¬ 
tering,  i.e.,  the  particles  in  a  Haze  M  distribution  may  just  as  well 
be  considered  as  Rayleigh  scatterers  for  wavelengths  greater  than  10p. 
It  is  only  for  the-  larger  particles,  cf.  figures  80  through  84,  which 
occur  with  small  probability,  that  we  see  any  appreciable  peakedness 
in  the  forward  direction.  The  cumulative  probability  function  also 
indicates  the  near  isotropy  at  \  =  10. Op.  Here  the  probability  of 
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Figure  55.  Haze  M  Distribution 
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SCATTERING  ANGLE  (DEG) 


Figure  56,  Particle  Scattering  Function  for  x  =  0.5u 
and  x  =  1.06,  Haze  M  Distribution 
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Figure  59.  Particle  Scattering  Function  for  X  =  0.5u 
and  x  =  16.13,  Haze  M  Distribution 
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Figure  61.  Volume  Scattering  Function  I j (©)  for  Haze  M 
Distribution  at  X  =  0.5u 
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Figure  63. 


Cumulative  Distribution  for  Haze  M 
Distribution  at  ).  -  d.5u 
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Figure  65.  Particle  Scattering  Function  for  \  =  2. On 
and  x  =  1.02,  Haze  M  Distribution 
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Figure  69.  Volume  Scattering  Function  Ii(e)  for  Haze  M 
Distribution  at  \  =  Op 
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SCATTERING  FUNCTION  I2(e) 


Figure  70.  Volume  Scattering  Function  I2(e)  for  Haze  M 
Distribution  at  X  =  2. Op 
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Cumulative  Distribution  for  Haze  M 
Distribution  at  \  =  2.0m 
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Figure  74.  Particle  Scattering  Function  for  X  =  6.0y 
and  x  s  0.76,  Haze  M  Distribution 


FiguK  75.  Particle  Scattering  Function  for  \  =  6. Op 
and  x  =  1.34,  Haze  M  Distribution 
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Figure  76.  Particle  Scattering  Function  for  \  =  6. Op 
and  x  =  2.09,  Maze  M  Distribution 
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Figure  79.  Cumulative  Distribution  for  Haze  M 
Distribution  at  x  =  6.0y 


Figure  80.  Particle  Scattering  Function  for  X  =  10. Op 
and  x  =  0.0G,  Haze  M  Distribution 
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Fiqure  83.  Particle  Scattering  Function  for  \  =  10. On 
and  x  =  0.81,  Haze  M  Distribution 
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Figure  84.  Particle  Scattering  Function  for  X  =  10. On 
and  x  =  1 .26,  Haze  M  Distribution 
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scattering  through  an  angle  between  0°  and  71°  is  the  same  as  that 
for  angles  between  0°  and  21°  for  the  cumulus  distribution  and 
A  =  10.0m.  One  can  see  from  figures  85  and  86  that  the  normalized 
scattering  function  at  A  =  10. On  rest.:, die  scattering  functions  for 
Rayleigh  scattering.  Figure  87  shows  the  almost  uniform  cumulative 
probability  functions.  Table  2  contains  the  extinction  coefficients 
for  the  Haze  M  distribution  as  defined  by  equation  (324).  We  note 
that  the  extinction  coefficient  decreases  with  an  increase  in  w,  'e- 
length,  i.e.,  red  is  attenuated  less  than  the  shorter  wavelengths  as 
was  first  shown  by  Rayleigh. 

Table  2 

EXTINCTION  COEFFICIENTS  FOR  THE  HAZE  M  DISTRIBUTION 


*(u) 

^ext 

0.5 

1.001  x  1.0“ 6 

1.0 

9.083  x  10“7 

2.0 

4.297  x  10“7 

4.0 

1.222  x  10"7 

6.0 

4.040  x  10“8 

8.0 

1.677  x  10“8 

10.0 

7.931  x  I 0“9 

4.  DEPOLARIZATION  BY  A  CLOUDY  ATMOSPHERE  AT  0.5,  4.0,  and  10. On 
In  this  subsection  we  consider  the  propagation  of  an  electromag¬ 
netic  wave  through  an  utmosphere  containing  cumulus  clouds.  The  photon 
source  will  be  pointed  in  the  downward  direction  (see  figure  88)  from 
a  fixed  altitude  of  13  km.  We  assume  that  the  source  is  linearly 
parallel  polarized  and  that  the  source  beamwidth  ’s  2°  x  2°.  The 
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Figure  88.  Cumulus  Cloud  Geometry 
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altitude  of  the  cloud  bottom  is  10  km  while  the  altitude  of  the  cloud1 
top  is  11  km.  We  first  assume  that  only  Rayleigh  scattering  exists 
above  and  below  the  cloud,  then  only  Mie  scattering,  and  finally  we 
assume  that  no  scattering  occurs  outside  the  cloud.  Inside  the  cloud 
we  consider  only  Mie  scattering.  We  determine  the  depolarization  of 
the  linearly  polarized  wave  at  wavelengths  of  0.5,  4.0,  and  iO.Ou. 

The  ground  is  assumed  to  be  a  diffuse  scatterer  so  that  the  incident 
wave  is  completely  depolarized,  i.e.,  equal  parallel  and  perpendicular 
components  are  scattered.  We  initially  assume  a  ground  albedo  of  1.0, 
i.e.,  all  incident  energy  is  scattered  back  into  the  atmosphere,  then 
we  assume  that  no  ground  scattering  occurs.  The  source  is  placed  at 
an  altitude  of  13  km.  Two  isotropic  detectors  are  placed*  at  altitudes 
of  9  km  (Dl)  and  13  km  (D2)  with  D2  displaced  horizontally  by  0:3m 
from  the  source. 

In  simulating  the  problem  described  above  using  the  Monte  Carlo 
technique  we  use  5,000  and  10,000  histories.  Experience  gained  in 
running  several  problems  has  lead  to  the  selection  of  10  collisions 
(see  figure  109).  Actually,  as  we  can  see  from  figure  109,  at  6 
collisions  the  received  intensity  has  dropped  by  3  orders  of  magnitude 
and  we  could  probably  carry  the  number  of  collisions  out  to  only  6 
with  acceptable  results. 

Fiqure  89  shows  the  assumed  stratification  for  the  atmosphere. 

The  optical  depth  increases  linearly  with  altitude  in  steps  as  follows 

|10"5h  0  <  h  <  10  km 

10"3h  -  9.9  10  km  <  h  <  11  km  (335) 

6.1  x  10” 7h  +1.1  11  km  <  h  <  500  km 
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F  ure  89.  Atmospheric  Stratification 


The  optical  depth  of  the  cloud  is  1.0  as  can  be  seen  from  figure  89. 
Normally  the  atmosphere  is  stratified  exponentially  [15]  but  for  small 
distances  as  we  have  assumed  for  this  problem  linear  stratification  is 
a  good  assumption  for  the  type  of  atmosphere  selected.  This  basic 
stratification  will  be  scaled  up  for  simulation  of  greater  densities. 

As  we  have  stated  the  polarization  assumed  for  this  problem  is 
parallel.  This  means  that  the  electric  field  vector,  according  to  the 
convention  adopted  in  section  III,  is  in  the  vertical  plane  t.irough  the 
propagation  vector  and  the  vertical  axis  referred  to  the  source.  Figure 
90  portrays  our  polarization  convention.  The  source  is  located  at  the 
origin.  In  this  example  a  photon  is  being  emitted  from  the  source  in 
the  direction  (e,$),  the  electric  field  vector  is  in  a  vertical 
plane  containing  the  prorogation  vector  k  and  the  vertical  axis. 

Figure  91  shows  the  results  of  the  Monte  Carlo  simulation  for  this 
problem.  The  ordinate  is  the  polarization  factor  and  is  defined  as 

D  -  jCi  (336) 

r  i 
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Figure  90.  Polarization  Reference  Plane  for  Emitted  Intensity 

where  I  and  I  are  the  perpendicular  and  parallel  intensity  components, 
respectively.  The  abscissa  is  the  cosine  of  the  polar  angle  for  each 
detector.  Recall  that  the  polar  angle  is  measured  from  the  line  join¬ 
ing  the  detector  and  the  source  with  0°  toward  the  source  and  180°  away 
from  the  source  behind  the  detector.  Therefore  in  figure  91  both  detec¬ 
tors  are  "looking"  directly  at  the  source  when  the  polar  cosine  is 
euqal  to  +1.0  and  directly  away  when  the  polar  cosine  is  -1.0.  Detector 
2  is  only  0.3  meters  from  the  source  and  the  magnitude  of  the  intensity 
from  all  directions  is  such  that  the  variance  of  the  samples  is  less  than 
two  orders  of  magnitude  from  the  mean.  Detector  1  is  4  kilometers  from 
the  source  and  due  to  the  ncrrowness  of  the  source  beam  (2°  x  2°)  the 
variation  of  the  intensity  and  therefore  also  of  the  polarization  fac¬ 
tor  is  large  between  values  of  the  polar  cosine  equal  to  -0.7  and  +0.7, 
i.e.,  the  data  are  meaningful  only  for  forward  scattering  angles  from 
0°  to  45°  and  for  back  scattering  angles  from  135°  to  180°. 
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In  figure  92  we  show  the  depolarization  factor  for  the  same  problem 
except  at  A  =  0.5u.  In  comparing  figures  91  and  92  we  can  see  that  for 
detector  1  the  depolarization  is  not  quite  as  great  as  it  is  for  detec¬ 
tor-2  at  0.5y.  The  intensity  scattered  directly  back  from  the  cloud 
and  the  ground  to  detector  2  remains  depolarized  at  0 . 5y  and  at  10. Ou. 

Figures  93  and  94  show  the  same  data  as  shown  in  figure  92  except 
expanded.  Figure  93  gives  the  polarization  factor  for  de sector  1  for 
polar  angles  from  0°  to  2°,  i.e.,  forward  scattering  through  the  cloud 
and  figure  94'  gives  the  polarization  factor  for  polar  angles  from  89° 
to  91°,  i.e.,  backscattering  from  the  cloud. 

Figure  91  shows  that  for  detector  1  located  below  the  cloud  (see 
figure  68)  that  the  forward  scattered  energy  (polar  cosine  equal  to 
+1.0)  is  only  partially  depolarized  and  at  that  point 

D  *  0.8  (337) 

while  for  backscattering  (polar  cosine  equal  to  -1.0)  from  the  ground 
with  an  albedo  equal  to  T.O  the  wave  is  completely  depolarized,  i.e., 

Ir  =  It  (338) 

or 

D  =  0  (339) 

This  is  as  expected  since  we  are  assuming  that  the  ground  is  a  diffuse 
scatterer  and  completely  depolarizes  the  incident  wave. 

For  detector  2  we  can  see,  if  we  ignore  the  small  variations,  that 
the  received  intensity  is  almost  completely  depolarized  for  all  direc¬ 
tions.  The  cloud  and  ground  are  in  the  direction  where  the  polar  cosine 
is  equal  to  0.  For  figure  91  and  subsequent  data  the  received  intensity 
is  integrated  over  all  azimuth  angles,  but  the  magnitude  of  the  inten¬ 
sity  received  from  angles  other  than  forward  and  backscatter  can  be 
considered  to  be  negligible;  for  example,  for  detector  2  the  intensity 
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Figure  93.  Forward  Scatter  Polarization  Factor  for  Cumulus 
Cloud  at  A  =  0.5g  with  Ground  Scatter 
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Figure  94.  Backscatter  Polarization  Factor  for  Cumulus  Cloud 
at  X  =  0.5y  with  Ground  Scatter 
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received  from  an  azimuthal  interval  of  45°  while  pointing  directly  at 
the  cloud  4s  four  orders  of  magnitude  greater  than  the  intensity 
received  from  an  interval  of  45°  but  pointed  180°  from  the  cloud,  i.e., 
in  the  positive  vertical-  direction. 

The  data  just  presented  in  figures  91,  92,  93,  and  94  include,  as 
pointed  out,  reflections  from  the  ground  as  well  as  from  the  cloud.  In 
order  to;  determine  the  effect  of  the  ground  on  the  depolarization  of  the 
incident  wave  the  ground  was  assumed  to  be  highly  absoroent  and  an  albedo 
of  0.001  was  selected  compared  to  an  albedo  of  1,0  for  the  previous  fig¬ 
ures.  Figures  95  and  96  show  the  results  of  this  modification  for  5,000 
histories  and  a  =  10. 0y.  In  figure  95  we  can  see  that  the  average, 
depolarization  factor  for  detector  1  is  approximately  0.9  compared  to 
0.8  for  forward  scatter  in  figure  91.  The  approximate  average  polariza¬ 
tion  factor  for  detector  2  for  intensity  backscattered  from  the  cloud 
in  figure  96  is  approximately  0.7  while  if  ground  reflections  are 
included  as  in  figure  91  the  backscattered  intensity  is  completely 
depolarized.  One  can  therefore  conclude  that  at  a  =  10. On,  for  the  prob¬ 
lem  defiled,  the  diffuse  scattering  ground  will  cause  a  major  portion 
of  the  depolarization  of  the  incident  wave.  It  should  be  pointed  out 
that  the  data  obtained  for  figures  95  and  96  was  obtained  while  allowing 
up  to  10:  photon  collisions.  A  separate  problem  was  run  using  10,000 
histories  while  allowing  only  1  collision;  the  scattered  intensity  was 
not  depolarized  at  all  which  agrees  with  theory.  Therefore  we  can  see 
that  for  '  fairly  thin  cloud  with  an  optical  depth  of  1.0,  depolariza¬ 
tion  due  to  multiple  scattering  is  significant. 

In  pursuing  this  problem  further  we  determined  the  degree  of 
depolarization  for  thicker  cumulus  clouds  at  A  *  0.5 n  and  4. On.  The 
wavelength  4. On  was  selected  to  coincide  with  the  mode  radius  for  the 
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95.  Forward  Scatter  Polarization  Factor  for  Cumulus  Cloud 
at  X  -  10. Op  with  No  Ground  Scatter 
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Figure  96.  Backscatter  Polarization  Factor  for  Cumulus  Cloud 
at  X  =  10. Op  with  no  Ground  Scatter 
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size  distribution  of  the  cumulus  cloud  (figure  Z2)\  The  optical  depth 
for  x  =  0.5m  was  selected  to  be  3.0  while  at  x  =  4.0m  we  selected,  t  = 
4.0.  The  values  of  the  optical  depth  at  0.5m  and  4.0m  were  ecu  ad 
according  to  the  ratio  of  the  individual  extinction  coefficient  which 
are 


°ext  (4-0>l) 
“ext 


4 

3 


(340) 


The  scaling  in  the  optical  depth  t  follows  from  equations  c  .73)  and 
(274)  where  we  see  that  the  volume,  scattering  functions  a  ;e  proportional 
to  x2.  Any  differences"  in  t  due'  to  frequency  change  for  a  particular 
size  distribution  must  therefore  be  scaled  according- to  t.  /e  change  in 
extinction  coefficient.  This  scaling  allows  us  to  maintain  the  same 
physical  atmospheric  model.  Again  we  assume  a. highly  absorbent  ground. 
Figures  97  and  98  show  the  backscattered  and  forward,  scattered  polariza¬ 
tion  factors  at  x  =  0.5y,  respectively.  Figures  99  and  100  show  the 
backscattered  and  forward  scattered  polarization  factors  at  X  =  4.0m, 
respectively.  We  can  see  that  at  X  =  0.5m 

D  «  -0.65  (341) 


while  lit  X  =  4.0m 


D  »  -0.45 


(342) 


We  therefore  conclude  that  for  the  cumulus  cloud  distribution  that  a 
greater  amount  of  depolarization  occurs  at  x  =  4.0m  than  at  x  =  0.5m. 
This  can  be  explained  by  noting  that  the  scattering  functions  for 
cumulus  cloud  distributions  are  not  as  highly  peaked  in  the  forward 
and  backward  directions  at  larger  wavelengths  (figures  36  and  37)  as 
they  are  at  the  smaller  wavelengths  (figures  28  and  29).  This  change 
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Figure  97.  Backscatter  Polarization  Factor  for  Cumulus  Cloud 
at  X  =  0.5y  with  no  Ground  Scatter 
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Figure  98.  Forward  Scatter  Polarization  Factor  for  Cumulus  Cloud 
at  X  =  0,5u  with  no  Ground  Scatter 
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Figure  99.  Backscatter  Polarization  Factor  fot 
at  A  =  4. On  with  no  Ground  Scattei 
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Figure  100.  Forward  Scatter  Polarization  Factor  for  Cumulus  Cloud 
at  X  =  4.0y  with  no  Ground  Scatter 


in  scattering  function  results  in  a  greater  degree  of  multiple  scat¬ 
tering  for  longer  wavelengths  and  thus  causes  a  smaller  polarization 
factor  at  that  wavelength  compared  to  that  of  0.5p. 

The  forward-scattered  polarization  factors  at  X  =  0.5p  and  x  =  4. Op 
are  shown  in  figures  98  and  100  and  we  can  see  that  the  average  value 
for  both  wavelengths  is  approximately 

D  «  -0.6  (343) 

and  we  therefore  conclude  that  the  forward  scattered  polarization  is 
not  as  sensitive  to  frequency  as  is  the  backscattered  intensity. 

We  continue  to  isolate  the  cause  of  depolarization  by  not  allowing 
either  ground  or  Rayleigh  scattering  at  X  =  4. Ops-  Figure  101  is  the 
polarization  factor  for  backscatter  from  the  cloud;  we  note  that  now 
D  ss  -0.65  versus  D  «  -0.58  in  -figure  99'  when  we  include  Rayleigh  scat¬ 
tering.  Figure  102  shows  the  polarization  factor  for  forward  scatter 
from  the  cloud  for  no  ground  or  Rayleigh  scattering.  Comparing  figure 
102  with  figure  100  we  note  that  at  0°  D  «  -0.88  in  figure  102  versus 
D  «  -0.68  in  figure  100.  We  now  conclude  that  Rayleigh  scattering  causes 
an  increase  in  depolarization.  We  should  note  that  although  we  have' 
eliminated  Rayleigh  scattering,  we  have  replaced  it  with  Mie  scattering 
particles  outside  the  cloud  taken  from  the  cumulus  distribution. 

Now  we  completely  eliminate  any  form  of  scattering  outside  the  cloud. 
Figure  103  shows  the  backscatter  polarization  factor  from  the  cloud  and 
figi^re  104  the  forward  scatter  polarization  factor.  On  comparing  fig¬ 
ures  101  and  102  with  figures  103  and  104  we  note  that  the  change  is 

v 

small. 

We  now  compare  the  polarization  factors  for  no  scattering  outside 
the  cloud  at  X  =  0.5p  and  X  =  4. Op.  Figure  105  shows  the  backscatter 
polarization  factor  at  X  =  0.5p  and  we  note  that  at  90°  D«  -0.78  versus 
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Figure  102.  Forward  Scatter  Polarization  Factor  for  Cumulus  Cloud 
at  X  =  4, Op  with  no  Ground  or  Rayleigh  Scattering 


Figure  103.  Backscatter  Polarization  Factor  for  Cumulus  Cloud 
at  X  =  4.0m  with  no  Scattering  Outside  Cloud 
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Figure  104.  Forward  Scatter  'olarization  Factor  for  C 
at  \  =  4 . Oja  with  no  Scattering  Outside  Cloud 


Figure  105.  Backscatter  Polarization  Factor  for  Cumulus  Cloud 
at  A  =  0.5y  with  no  Scattering  Outside  Cloud 


D  «  -0.65  at  A  -  4. Op  in  figure  101.  Figure  106  shows  the  forward  scat¬ 
ter  polarization  factor  at  A  =  0.5p  and  at  0°  D  «  -0.8  versus  D  «  t.0.78 
at  A  =  4. Op  in  figure  102.  We  therefore  conclude  that  the  degree  of 
depolarization  is  greater  at  A  =  4. Op  than  -at  A  =  0.5p  for  a  Cumulus 
Cloud  distribution.  Ther  percentage  change  in  polarization  factor  from 
A  =  0.5p  to  A  =  4. Op  is  approximately  9  percent  for  a  cloud  of  medium 
thickness  (t  =  4.0). 

5.  DEPOLARIZATION  BY  A  HAZY-CLOUDY  ATMOSPHERE  AT  0.866p 

Consider  now.  the  geometry  shown  in  figure  107.  Again  as  in  the  last 
problem  we  have  a  source  and  two  detectors.  Here  detector  1  is  located 
just  10  meters  off  the  ground.  The  source  and  detector  2.  are  located 
at  an  altitude  of  5  km  with  detector  2  located  0.5  meter  horizontally 
from  the  source.'  The  source  teamw.idth  is  2°  x  2°  and  emits  photons  at 
A  =  0.866p.  The  atmosphere  is  composed  of >Haze  M  with  a  cumulus  cloud 
between  1  km  and  3  km  from  the  ground.  For  this  problem  the  ground 
albedo  is  assumed  to  be  0.8. 

Here  we  assumed  that  the  atmosphere  is  exponentially  stratified  and 
the  optical  depth  is  given  by  (see  section  V)  [24] 


x(h)  =  0O  h(i  -  e‘h/H)  (344) 

where 

aQM  =  1.02  x  10"6  for  Haze  M 
°0C  =  1*70  x  10-4  for  cumulus  clouds 
Hm  =  0.98  km  for  Haze  M,  and 
Hc  =  2.00  km  for  cumulus  clouds 
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Jhe  value  H--  =  2.00  km  was  assumed  while  oQM,  aQC,  and  are  measured 
values  [15]..  Figure  108  is  a  plot  of  equation  (344)  for  the  .geometry 
shown  in  figure'107.  For  this  model  the  optical  depth  of  the  atmosphere 
is  -1.7.  .The  Size -particle  distributions  for  Cumulus  Clouds  and  Haze  M 
are  given  by  equations  (311)  and  (325),  respectively.  We  assume  that 
Rayleigh  scattering  is  negligible  everywhere  so  that  only  Mie  scatter¬ 
ing  occurs.  Up  to .three  reflections  from  the  ground  were  allowed. 

In  figure  109  we  have  shown  the  relative  intensity  at  detector  1 
as  a  function  of  the  number  of  collisions  during  a  single  history  for 
10,000  histories.  In  the  simulation  we  consider  up  to  10  collisions 
and' we  can  see  thit  the  intensity  becomes  negligible  for  collisions 
greater  than,  say,  5  collisions. 

Figures  110  and  111  represent  the  depolarization  factor  as  a  func¬ 
tion  of  the  polar  cosine  where  the  polar  angle  is  defined  as  in  sub¬ 
section  4.  From  figure  110  we  can  see  that  for  detector  1  which  is 
located  directly  below  the  source  and  10  meters  above  the  ground  (see 
figure  107)  the  received  intensity  is  only  slightly  depolarized  when 
"looking"  directly  at  the  source,  i.e., 

0  =  -0.8  (345) 

for  polar  cosine  equal  to  +1.0.  On  the  other  hand  the  intensity  is 
completely  depolarized  for  scattering  from  the  ground  at  a  polar 
cosine  of  -1.0.  Depolarization  increases  gradually  reaching  a  value 
of  zero  at  approximately  90°  or  polar  cosine  equal  to  0. 

For  detector  2  which  is  located  next  to  the  source  the  received 
intensity  is  almost  completely  depolarized,  i.e.,  D  «  0.  The  problem 
was  then  rerun  with  no  ground  reflections  to  influence  the  polarization 
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Figure  110.  Polarization  Factor  for  Cumulus/Haze  M  for  D1 
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Figure  m.  Polarization  Factor  for  Cumulus/Haze  M  for  D2 
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factor..  For  this  case  +he  backscattered  intensity  from  the  cloud  at  D2 
was  found  to  remain  parallel  polarized,  i.e., 

D  =  -1  (346) 

6.  DEPOLARIZATION  WITHIN  A  CUMULUS  CLOUD  AT  0.5:  and  10. On 

Let  us  now  investigate  the  problem  of  depolarization  of  a  parallel 
polarized  wave  when  both  the  source  and  detectors,  are  within  a  cumulus 
cloud  (figure  112).  The  cumulus  cloud  size  distribution  as  given  by 
equation  (311).  Here  again  we  will  assume  that  the  source  emits  photons 
in  a  2°  x  2°  beam  which  is  pointed' horizontally  in  line  with  two  iso¬ 
tropic  detectors.  The  source  and  detectors  are  located  at  an  altitude 
of  10  km.  The  cloud  extends  from  the  ground  up  to  a*-  altitude  of  20  km 
above  which  no  scatterers  exist.  The  ground  albedo  was  selected  to  be 
small  (0.1)  in  order  to  simulate  a  highly  absorbent  ground.  We  again 
allowed  up  to  three  reflections  from  the  ground.  Detector  D1  is  0,3 
meters  from  the  source  while  detector  D2  is  1  km  from  the  source.  With 
D1  we  will  measure  the  backscattered  perpendicular  and  parallel  polariza 
ti oh  components  while  with  D2  we  will  measure  the  perpendicular  and 
parallel  forward,  scattered  comf-.ients.  The  optical  thickness  of  the 
cloud  is  assumed  to  be  uniformly  increasing  with  altitude  as  shown  in 
figure  113.  This  variation  of  optical  depth  with  altitude  can  be 
represented  by  equation 

x(h)  =  5  x  10"4  h  (347) 

Figure  114  shows  the  result  at  10. On;  we  can  see  that  the  back- 
scattered  intensity  at  detectc*'  D1  is  slightly  more  than  half  depolar¬ 
ized  and  we'hawe 

D  =  -0.4  (348) 
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Geometry  for  Transmit-Receive  System  within  a  Cumulus  Cloud 


At  detector  D2  the  received  intensity  is  completely  depolarized. 

The  complete  depolarization  at  D2  was  later  found  to  be  caused  by 
ground  scatter.  We  can  see  that  for  D2  the  polarization  factor  begins 
to  vary  considerably  for  polar  cosines  less  than  0.9  («26°).  Here, 
again,  we  have-allowed  10  collisions  per  history  and  the  photon  inten¬ 
sity  is  attenuated  quite  rapidly  as  it  propagates  toward  D2  from  the 
source.  Rayleigh  scattering  was  not  allowed  for  this  problem. 

Figure  115  shows  the  polarization  factor  for  detector  1  for  the 
backscattered  intensity  at  10. Oy.  Notice  that  the  polar  angle  covers 
the  interval  from  177°  to  180°,  therefore  for  direct  and  quasi -direct 
backscatter  the  intensity  is  not  completely  depolarized  and  agrees  with 
the  larger  angle  interval  shown  in  figure  114. 

Figure  116  is  an  expanded  version  of  figure  114  for  detector  2  in 
the  Interval  from  0°  to  2°;  here  one  can  see,  if  one  neglects  the  small 
statistical  fluctuations,  that  the  wave  is  completely  depolarized. 

To  determine  the  degree  of  depolarization  caused  by  ground  reflec¬ 
tions  and  scattering  within  the  cloud,  we  eliminated  all  ground  reflec¬ 
tions  and  found  that  for  t  =  10.0  the  received  intensity  was,  for  all 
practical  purposes,  still  parallel  polarized,  i.e.,  D  =  -1  for  parallel 
polarization.  We  thus  conclude  that,  for  this  specific  problem,  the 
depolarization  is  caused  by  ground  scatter. 
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Figure  116.  Forward  Scatter  Polarization  Factor  for  Cumulus 
Cloud  at  \  =  10. Ou  and  t  =  10.0 


SECTION  VI 

DISCUSSION  AND  CONCLUSIONS 

In  this  study  we  have  set  out  to  determine  the  effect  of  scattering 
on  waves  propagating  through  the  atmosphere.  We  have  centered  our  atten¬ 
tion  in  the  infrared  part  of  the  spectrum  although  we  have  compared  our 
results  to  A  =  0.5y  which  ts  in  the  visible  portion  of  the  spectrum. 

In  section  I,  we  presented  a  brief  discussion  of  the  various  mech¬ 
anisms  which  may  cause  the  depolarization  of  electromagnetic  waves 
including  nonlinear  effects.  We  stated  that  although  polarization  has 
been  used  in  extracting  information  concerning  scattering  tragets,  prob¬ 
lems  in  controlling  the  polarization  state  of  an  electromagnetic  wave 
are  usually  encountered  especially  at  the  microwave  frequencies.  We 
also  pointed  out  that  in  separating  the  intensity  scattered  back  from 
multiple  targets  the  effects  of  multiple  scattering  must  be  taken  into 
account  for  dense  media. 

In  sections  II  and  III  we  derived  the  transfer  equation  for  polar¬ 
ized  waves  and  showed  how  one  might  treat  a  polarized  wave  of  arbitrary 
polari.  .tion.  In  section  II  we  gave  an  indication  on  how  we  would 
approach  the  modeling  of  scattering  of  polarized  waves  in  the  atmosphere. 
In  section  III  we  derived  the  scattering  matrixes  for  Mie  and  Rayleigh 
scattering  which  we  subsequently  used  in  modeling  the  scattering  pro¬ 
cess  in  the  atmosphere.  All  scattering  particles  in  the  atmosphere  were 
assumed  to  be  either  Rayleigh  or  Mie  particles,  Rayleigh  particles 
being  those  less  than  a  wavelength  in  size  while  Mie  particles  being 
spheres  of  arbitrary  size. 
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In  section  IV  we  presented  the  simulation  model.  We  derived  the 
probability  functions  which  characterize  the  path  of  a  photon  in  the 
atmosphere.  We  showed  that  through  the  Monte  Carlo  simulation  of  a 
large  number  of  photon  histories  we  could  obtain  approximate  numerical 
results  to  the  scattering  of  electromagnetic  waves  in  the  atmosphere 
including  the  effects  of  multiple  scattering. 

In  section  V  we  obtained  some  results  for  specific  problems.  We 
used  the  Cumulus  Cloud  and  Haze  M  size  distributions  to  obtain  the 
volume  scattering  functions  for  incident  perpendicular  and' parallel 
polarization.  Using  the  volume  scattering  functions  we  modeled  cloudy 
and  hazy  atmospheres  for  different  optical  depths  and  wavelengths.  In 
subsection  4  of  section  V  we  saw  that  under  the  assumption  of  a  diffuse 
scattering  ground  which  caused  complete  depolarization,  li.e.,  equal 
parallel  and  perpendicular  components,  \the  backscattered  intensity  was 
completely  depolarized  while  the  forward  scattered  intensity  was  only 
slightly  depolarized.  We  showed  this  to  be  true  at  wavelengths  of  0.5, 
4.0,  and  IO.Ojj.  We  then  saw  that  if  we  considered  only  the  intensity 
scattered  from  the  cloud  that  there  was  no  depolarization  for  single 
scattering  and  a  significant  amount  of  depolarization  for  multiple 
scattering.  By  comparing  the  backscattered  polarization  factor  at 
wavelengths  of  0.5  and  4.0y  for  the  same  physical  cloud,  we  saw  that 
the  depolarization  was  greater  at  4.0y  than  at  0.5y.  A  comparison  of 
the  volume  scattering  functions  at  these  two  frequencies  (figur*5*  .'8 
and  44)  reveals  that  this  result  is  to  be  expected  since  the  fu  {-tion 
width  increases  with  wavelength. 

In  subsection  5  of  section  V  we  modeled  a  low  lyir,.  cumulus  cloud 
in  a  hazy  atmosphere.  The  depolarization  for  this  model  is  comparable 
to  that  discussed  above.  In  addition  we  showed  how  the  relative 
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scattered  intensity  varies  as  a  function  of  collision  number.  The 
results  shown  in  figure  109  compare  favorably  to  those  obtained  by 
Bruning  and  Lo  [29]  for  multiple  scattering  of  electromagnetic  waves 
by  spheres. 

Thus,  based  on  the  results  of  this  study  we  are  lead  to  the  con¬ 
clusion  that  the  degree  of  depolarization  of  electromagnetic  waves  by 
scattering  In  the  atmosphere  can  be  obtained  through  a  Monte  Carlo 
simulation  similar  to  the  one  used  in  this  study.  We  have  shown  that 
the  polarization  state  of  an  electromagnetic  wave  could  be  used  to 
discriminate  between  a  diffuse  scattering  ground  and  a  cloud  or  haze 
if  the  cloud  or  haze  is  not  "too  dense."  By  "too  dense"  we  mean  - 
optical  distances  of  approximately  3.0  per  km  or  greater,  i.e.,  it 
appears  that  these  3  collisions  per  kilometer  are  sufficient  to  cause 
significant  depolarization  of  the  incident  wave  for  the  wavelengths 
considered  here. 

Much  work  remains  to  be  accomplished  in> order  t'1  fully  characterize 
the  scattering  phenomena  in  the  atmosphere.  Measurements  of  polariza¬ 
tion  factors  at  various  frequencies  under  controlled  conditions  must 
be  made  to  verify  the  results  presented  in  this  paper.  The  value  of 
a  computer  simulation  lies  in  the  ability  in  obtaining  accurate  results 
under  prescribed  conditions.  Once  a  simulation  model  provides  accurate 
results  the  economic  benefits  through  the  reduction  of  the  required 
experimental  measurements  are  obvious. 
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APPENDIX 


DERIVATION  OF  THE  TRANSFORMATION  MATRIX 


The  purpose  of  this  appendix  is  to  derive  the  transformation  matrixes 

»v 

of  F  end  F.  Here,  we  will  use  the  same  notation  used  by  Van  De  Hulst  in 
his  excellent  book  [14]. 

Let  the  transformation  matrix  F  be  defined  by 
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=  F 
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(1) 


Where  (I,  Q,  U,.  V)  are  the  scattered  Stokes  parameters,  and  Qq, 
Uq,  Vq)  are  the  incident  Stokes  parameters. 

In  general  the  scattered  electric  field  components  are  given  by 
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where  the  A^  constitute  a  transformation  matrix  of  (e|,  E^)  into  (e^, 
Er> 

The  Stokes  parameters  are  defined  by 
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Using  (2)  we  can  find  each  term  of  (-3)  as  follows: 
1  =  (A2  EI  +  A3  4  +  A3  Er) 
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where  the  last  two  terms  within  each  bracket  were  added  and  subtracted. 
We  can  new  write 


5  VZ  |a,I2  +  |A?|2  +  |A3P  +  |A4|2]  ;[Ef  EJ*  +  E l  E**] 

+  V2  i  |A,I2  +  lA2l2  -  lA3l2  +  M2]  K  Ef  -  £'r  Ef] 

+  W2  bl  \  +  A>„  +  A2  A*  +  A*  A3]  [E1  ?  +  Er] 


*  1/2  I:  A‘  **  +  A‘  +  *2  At  -  A*  *>]'  [E1  E’*  -  C  Erj 


or  (cf.  3) 
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Using  the  same  procedure  for  Qs  U,  and  V,  we  obtain 

Q=  1/2  [-  |AJ2  +  | A2 1 2  +  |A3 1 2  -  |AJ2]I0^ 
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+  1/2  [-  A1  A*  -  A*  A4  +  A2  A*  +  A*  A3]  UQ 

-  1/2  ^  aJ  -  A?  At*  +  A2  A*  -  A2  A,]  Vq 


(5) 
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1/2  |Aj  A3  +  Aj  A3  +  A2  A4  +  A2  A4J  Iq 


*  -* 


+  1/2  [-  Aj  A3  -  Aj  A3  +  A2  A4  +  A2  A4J  Qq 

+  i/2  [a,  a;  *  A*  a2  +  a3  a*  +  A^  a  J  u0 

-  1/2  [-  Aj  A2  +  Aj  A2  -  Aj  A4  +  A3  A4  j  Vq 


A*  +  A*  A2  +  A3  -A*  +  A*  A4]  U( 


V  =  1/2  [-  Aj  A3  +  Aj  A3  +  A2  A4  -  A2  A4j  1^ 
+  1/2  [Aj  A*  -  A*  A3  +  A2  A*  -  A*  A4  j Qq 
+  1/2  Aj  A2  +  Aj  A2  +  A3  A4  -  A3  A4]  Uj. 

+  "2  [A,  <  +  <  *.  -  A3  <  -  "a  \]  V0 


*  * 


We  now  make  the  following  definitions 

Mk  =  lAl|f 

M’ 

*»  ■  SJ1/2  (Ai  <  +  Ai  m 

■ -Dki  ^  (aj  <  -  ^  Ak) 


for  j  =  1,  2,  3,  4  and*k  =  1®»  3,  4.  Substituting  (8)  into  (4) 
through  (7)  results  in  M 


yT&B&Sfi"*  Jim  j 


while 


F  = 


/ 


M2  ,  Mg  , 

M4  ,  Mj  , 


S 

°23  * 

5  >1 » 


^23\ 


-  D 
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2S24j  2S31 ,  S21  +  S34,  -  D21  +  D34 
^2D24,  2D31,  D21  +  D34,  S21  -  S34  J 

Equations  (9)  and  (10)  are  the  most  general  transformations.  For 
spherical  particles  A3  =  A4  =  0;  then  (?)  and  (10)  simplify  to 


(10) 


F  = 


M1 

+  M2 

)-1/2( 

“Mi 

),  0  , 

0  \ 

+  m2 

)•’«( 

Mi 

+  M  ' 

j,  0  , 

o  1 

0 

> 

0 

»  S21*  ■ 

0 

0 

*  °21» 

S2./ 

(ID 


and 


M2,  0  ,  0  ,  0  \ 

0,MP0  ,  0 

0  »  0  ,  S21, :  -  D2j 

\  o  >  o  i  d21  ,  s21  ji 


(12) 


Equations  (11)  and  (12)  are  the  desired  transformation  matrixes. 


\ 
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